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ABSTRACT 


A  COMPUTER  MODEL  FOR  ANALYSIS  OF  DIELECTRICALLY 
ENHANCED  SINGLE-ZONE  TUNNEL  DRYERS 

by 

WILLIAM  MELENDEZ,  B.S. 

SUPERVISING  PROFESSOR:  PHILIP  S.  SCHMIDT,  Ph.D. 


Dielectric  heating  was  incorporated  into  the  existing  continuous 
single-zone  convective  dryer  model.  The  energy  balance  equation  for  a  control 

volume  around  the  drying  product  was  nondimensionalized,  and  the  resulting 

9 -  - 

dimensionless  parameters  were  used  to  generate  a  unique  characteristic  drying 
curve  and  to  define  drying  regimes  for  dielectrically  enhanced  drying.  Numerical 
models  were  developed  to  represent  the  heat  generation  inside  the  wet  product  and 
the  position  of  the  evaporative  front  as  a  function  of  moisture  content  The 
governing  equations  were  modified  to  account  for  the  thermal  influence  of  vapor 
flux  in  the  dry  region.  The  enhanced  computer  code  is  capable  of  predicting  the 
drying  rate  and  temperature  history  for  the  product  and  calculating  the  required 
dryer  length,  and  the  drying  time  for  dielectrically  enhanced  drying.  ^ u 
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CHAPTER  1 


INTRODUCTION 

The  use  of  dielectric  heating  in  drying  processes  has  been  an  area  of 
research  interest  for  several  years.  A  number  of  papers  have  been  published  on  the 
enhancing  effects  of  dielectric  drying.  Research  and  industrial  experience  has 
shown  that  dielectrically  enhanced  drying  can  reduce  thermal  degradation, 
producing  a  better  quality  product  than  conventional  drying  methods.  In  addition, 
dielectric  heating  allows  the  use  of  smaller  ovens,  attains  faster  drying  rates,  and 
maintains  a  cleaner  environment  relative  to  conventional  drying.  The  combination 
of  these  factors  can  result  in  a  more  efficient  drying  plant. 

Despite  the  high  efficiency  and  operational  advantages  of  dielectric 
drying,  its  adoption  in  industry  has  been  slow  for  two  main  reasons:  1)  high  capital 
cost  of  dielectric  heating  equipment,  and  2)  lack  of  predictive  tools  to  assess  the 
best  dielectrically  assisted  drying  configuration  for  a  particular  product 
Incorporation  of  dielectric  heating  involves  a  relatively  large  capital  investment,  and 
therefore,  must  be  judiciously  evaluated  and  optimized  for  its  successful 
implementation. 

1.1.  DRYER  MODELS 

In  the  pursuit  of  more  efficient  dryer  operation  and  less  expensive  dryer 
development,  process  engineers  have  developed  models  with  the  objective  of 
predicting  the  behavior  of  a  particular  dryer  under  specified  operating  conditions. 
As  documented  by  Reay  [1],  there  are  three  types  of  dryer  models  :  1)  the  scale-up 
factor  model  used  to  design  dryers;  2)  the  input/output  model  used  to  operate 
dryers;  and  3)  the  fundamental  model  based  on  the  fundamental  physics  of  the 
drying  process.  The  scale-up  factor  model  is  used  when  pilot  plant  tests  are 
conducted  and  the  operating  conditions  are  varied  until  product  specification  is  met; 
then,  the  dimensions  of  the  full-scale  dryer  are  estimated  using  empirical  scale-up 
rules  based  on  past  experience.  The  input/output  model  is  based  on  an  intuitive 
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understanding  gained  from  learning  the  impact  of  changes  on  product  and  operating 
conditions  through  a  posteriori  correlation  of  inputs  and  outputs.  The  fundamental 
model  is  based  on  a  mathematical  representation  of  the  physical  processes  occurring 
in  the  dryer,  but,  it  may  also  require  some  experimental  data  on  the  characteristics 
of  the  product.  The  scale-up  factor  and  input/output  models  cannot  be  extrapolated 
outside  the  range  of  past  experience;  the  fundamental  model,  however,  should  be  a 
more  powerful  tool  since  its  range  of  application  is  not  circumscribed  by  past 
experience  [  1  ]. 

The  dryer  model  discussed  in  this  report  is  of  the  fundamental  type.  A 
fundamental  dryer  model  must  have  two  components:  a  product  model  and  an 
equipment  model.  Figure  1.1  depicts  the  relationship  between  these  two  model 
elements.  The  product  model  characterizes  the  drying  kinetics  of  the  product 
independently  of  the  dryer  type.  Its  parameters  can  be  attained  from  bench  scale 
tests  performed  on  a  small  quantity  of  product  The  function  of  the  product  model 
is  to  determine  drying  rates  at  any  moisture  content,  temperature,  dielectric  power, 
or  local  drying  air  conditions.  Three  elements  are  necessary  for  a  complete  product 
model:  1)  a  drying  rate  function  for  a  small  quantity  of  the  product  as  a  function  of 
temperature,  air  psychrometric  properties  and  dielectric  heating  rate;  2)  an  algorithm 
for  predicting  the  drying  rate  at  any  given  operating  condition  and  moisture  content; 
3)  information  on  hygroscopicity  and  thermophysical  changes  of  the  product  with 
respect  to  its  moisture  content. 

The  equipment  model  represents  the  interaction  between  the  product 
model  and  any  energy  source  (air,  dielectric,  hot  plate)  independently  of  the 
particular  product.  This  model  relies  on  the  product  model  for  the  numerical  values 
of  drying  rate  and  product  properties.  The  equipment  model  must  be  able  to 
determine  the  heat  transfer  rate  from  the  air  to  the  product,  calculate  the  change  in 
product  and  air  moisture  content  and  the  necessary  dryer  length  or  drying  time 
needed  to  accomplish  the  desired  final  product  moisture  content. 


Product  Model 


Drying  Rate  \ 

Product: 

Thermophysical  Propertie 
Size 

Density  J 


Moisture  Content  > 
Product  Temperature 
Air  Temperature 
Air  Humidity  > 


Figure  1.1  Components  and  elements  of  a  fundamental  dryer  model  [1]. 
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Figure  1.2  Continuous  convective  tunnel  dryer.  Adapted  from  Reference  [2]. 
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The  University  of  Texas  dryer  model  and  its  computer  code  algorithm  are 
the  results  of  efforts  to  develop  a  fundamental  model  of  the  drying  process.  This 
model  was  documented  by  Rake  [2]. 

The  product  model  is  based  on  the  characteristic  drying  curve  principle,  as 
predicated  by  Keey  [3,4, 5,6]  for  convective  drying.  The  equipment  model 
interfaces  the  drying  process  with  the  product's  drying  kinetic  model  in  a 
continuous  convective,  single- zone  tunnel  dryer,  which  is  shown  in  Figure  1.2. 

The  existing  dryer  model  is  able  to  handle  drying  with  concurrent  or  countercurrent 
air  flow  and  with  isothermal  or  adiabatic  air  conditions.  The  convective  heat 
transfer  and  mass  transfer  coefficients  are  estimated  from  Nusselt  number 
correlations.  The  equipment  model  predicts  all  the  air  conditions  in  the  dryer, 
including  the  wet  bulb  temperature,  adiabatic  saturation  temperature  and  regions 
where  condensation  is  imminent 

There  are  certain  drying  conditions  that  Flake's  dryer  model  cannot 
handle.  First,  the  code  cannot  handle  boiling  point  drying  (i.e.,  conditions  where 
the  product  reaches  saturation  temperature).  Second,  though  provisions  were  made 
to  account  for  internal  heat  generation,  a  product  model  with  dielectric  heating  is  not 
available.  The  effect  of  dielectric  heating  on  the  characteristic  drying  curve  is 
unknown  since  Keey  formulated  the  characteristic  drying  curve  for  convective 
drying  only.  These  limitations  form  the  basis  for  the  objectives  of  this  research. 

1 .2.  OBJECTIVES  OF  THE  PRESENT  WORK 

The  major  long-range  objective  of  the  Dielectric  Heating  Research 
Program  at  the  University  of  Texas  at  Austin  is  to  develop  a  complete  predictive 
capability  for  dielectric  heating  applications  in  industrial  heating  and  drying 
processes.  One  step  toward  this  goal  was  the  development  of  a  fundamental  dryer 
model  for  a  continuous  convective  tunnel  dryer,  and  a  computer  code  capable  of 
interfacing  the  tunnel  dryer  model  with  the  product  kinetics  model  [2].  The  overall 


objective  of  the  present  study  is  to  incorporate  dielectric  heating  into  the  existing 
dryer  model.  To  achieve  this  goal,  the  effects  of  dielectric  heating  on  the  existing 
product  model  and  equipment  model  need  to  be  explored  in  detail. 


The  specific  objectives  of  the  present  work  are: 

•  Determine  pertinent  parameters  to  normalize  dielectrically  assisted  drying  curves 
and  produce  a  “characteristic  drying  curve”  which  accounts  for  the  effects  of 
dielectric  heating. 

•  Incorporate  the  characteristic  drying  curve  into  the  product  model  of  the  existing 
computer  code. 

•  Incorporate  the  heat  generation  term  into  the  equipment  model  of  the  existing 
computer  code. 

•  Demonstrate  the  application  of  the  proposed  normalization  scheme  to  a  limited  set 
of  experimental  data  for  dielectrically  enhanced  drying. 


Modify  the  existing  single-zone  computer  code  to  handle  dielectrically  assisted 
drying. 


CHAPTER  2 


THEORETICAL  BACKGROUND 


2.1  DRYING  RATE.  VAPOR  FLUX  EQUATION 


Keey  [3]  develops  a  detailed  derivation  of  the  drying  rate  from  a  wet 
surface  based  on  Fick's  law.  The  drying  rate  is  expressed  as 
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Equation  2. 1  gives  the  maximum  drying  rate  when  the  surface  humidity  is 
evaluated  at  the  wet  bulb  temperature  inside  the  dryer,  since  at  this  condition  Ys 
becomes  the  saturation  humidity.  Three  major  assumptions  are  necessary  to  derive 
Equation  2.1.  First,  the  density  of  the  vapor  and  the  diffusion  coefficient  are 
assumed  constant.  Second,  the  vertical  velocity  at  the  concentration  boundary  layer 
is  assumed  to  be  negligible.  Third,  the  difference  between  the  self-diffusivity  (Daa) 
and  the  limiting  interdiffusivity  (Dab)  is  assumed  to  be  negligible  [22].  In  the  case 
of  drying  with  dielectric  heating,  the  assumptions  inherent  in  Equation  2.1  are  no 
longer  valid.  Drying  rates  are  much  higher  in  dielectric  heating  due  to  higher  heat 
input;  hence,  the  convective  vapor  flux  will  affect  the  concentration  boundary  layer 
on  top  of  the  product  surface.  In  addition,  dielectric  drying  can  occur  at 
temperatures  well  above  the  wet  bulb  temperature,  and  in  fact  above  the  dry  bulb 
temperature  of  the  airstream.  When  Ta  >  100°C,  the  difference  between  the  self- 
diffusivity  and  the  limiting  interdiffusivity  is  no  longer  negligible  [22]. 

When  dielectric  heating  is  employed,  the  internal  drying  mechanisms  can 
become  more  complicated  since  other  vapor  transport  drivers,  besides  Fickian 
diffusion,  begin  to  dominate.  Pressure  gradients  and  thermal  gradients  may  play 
significant  roles  in  the  determination  of  the  total  moisture  transport  within  the 


The  characteristic  drying  curve  principle  was  first  introduced  by  van  Meel 
in  1958  [21].  Since  then,  its  application  to  industrial  drying  has  been  extensively 
studied  and  elaborated  on  by  many  researchers,  especially  R.B.  Keey  whose  work 
forms  the  basis  for  the  product  model  used  in  the  University  of  Texas  diyer  model. 
The  principle  of  the  characteristic  drying  curve  states  that,  for  a  given  material,  it  is 
possible  to  develop  a  normalized  drying  rate  curve  which  is  nearly  independent  of 
external  drying  conditions  [4]. 

2,2.2  Background 

The  theory  behind  the  characteristic  drying  curve  is  well-documented  by 
Keey  in  several  articles  [3,5,6];  a  review  of  key  concepts  follows.  To  derive  the 
drying  rate  curve  for  a  given  material  using  experimental  data,  the  slope  of  the 
mass-time  curve  obtained  at  fixed  conditions  of  air  temperature,  velocity  and 
humidity  is  calculated  and  plotted  against  the  moisture  content  as  presented  in 
Figure  2.1.  This  plot  of  Nv  vs.  X,  called  the  drying  rate  curve,  is  used  extensively 
in  the  design  and  evaluation  of  conventional  drying  equipment  [2]. 


In  many  drying  processes,  an  initial  constant  rate  period  occurs.  This 
constant  drying  rate,  NVo,  is  the  maximum  drying  rate  at  which  a  wet  surface 
reaches  the  wet  bulb  temperature  under  convective  conditions.  This  rate  is  used  to 
normalize  the  drying  rate  curve.  The  constant  drying  rate  period  is  represented  by 
the  wet  surface  model  which  assumes  that  evaporation  takes  place  only  at  the 
surface  of  the  product  and  that  surface  liquid  is  replenished  with  liquid  from  the 
interior  of  the  product  by  capillary  action. 


At  some  point  in  the  convective  drying  process,  capillary  action  is  not 
sufficient  to  keep  the  surface  wet,  and  dry  spots  appear  on  the  surface.  With  time, 
a  dry  region  develops  between  the  interior  of  the  wet  product  and  the  surface.  This 
dry  region  acts  as  a  resistance  to  vapor  transport,  and  the  drying  rate  decreases. 
During  this  period,  called  the  falling  rate  period,  the  evaporation  rate  falls  as  the 
moisture  content  decreases.  This  period  is  often  characterized  by  a  “receding 
moisture  front  ”  in  which  evaporation  occurs  at  the  surface  of  a  shrinking  wet  core 
topped  by  a  dry  shell  of  material  [1].  Some  materials  may  experience  more  than 
one  falling  rate  period  due  to  different  internal  moisture  migration  mechanisms. 

The  critical  moisture  content,  Xc,  is  defined  as  the  average  moisture 
content  of  the  product  at  the  transition  point  between  the  constant  drying  rate  and 
the  falling  rate  period.  As  stated  by  Keey,  Xc  is  not  a  physical  property  of  the 
material,  since  its  value  depends  on  the  thickness,  initial  drying  rate,  and  diffusivity 
of  moisture  inside  the  material.  However,  it  is  an  identifiable  feature  of  the  drying 
rate  curve  for  a  given  material  dried  under  fixed  conditions.  Xc  is  the  scaling  factor 
for  the  normalized  moisture  content,  <l>. 

The  drying  intensity,  Q,  is  defined  as  the  ratio  of  the  drying  rate  at  the 
surface  of  the  product  to  the  moisture  diffusion  rate  inside  the  material.  Its  value 
will  determine  if  a  receding  front  will  be  observed  during  the  drying  process.  If  Q 
is  much  greater  than  2,  a  sharp,  well-defined  front  will  sweep  down  the  product;  if 
Q  is  much  less  than  2,  a  sharp  front  will  not  be  observed  [3]. 

The  characteristic  drying  curve  is  defined  as  the  normalized  drying  rate, 

/r,  plotted  against  O,  as  indicated  in  Figure  2.1.  When  using  this  curve,  it  is 
assumed  that  a  change  in  external  conditions  will  have  the  same  proportional  effect 
on  the  falling  rate  period  as  on  the  constant  rate  period.  Keey  states  that  there  is  no 
single  characteristic  drying  curve  for  a  given  material  and  that  a  single  curve  will  be 
observed  only  when  Q,  the  heat  transfer  Biot  number,  Bi,  and  the  evaporative 
resistance  coefficient,  y,  assume  some  value  [6,16]. 
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Dielectric  heating  has  been  proposed  for  drying  applications  due  to  the 
high  affinity  of  water  for  electromagnetic  energy  at  microwave  and  radio 
frequencies.  Dielectric  heating  can  be  advantageous  under  the  following  conditions 
[11]: 


There  is  inadequate  heat  transfer  due  to  internal  resistance  to  heat  flow. 


There  is  inadequate  heat  transfer  due  to  a  small  external  heat  transfer  coefficient 


The  internal  mass  transfer  is  slow  due  to  the  moisture  being  held  in  the  pendular 
state  and/or  the  moisture  movement  is  limited  by  vapor  diffusion. 


The  resulting  quality  of  the  dried  product  is  poor  due  to  the  conventional  drying 
system. 


Dielectric  heating  serves  to  volumetrically  heat  a  material  by 
electromagnetic  energy  generated  at  radio  (3-300  MHz)  or  microwave  (300-3000 
MHz)  frequencies.  This  form  of  heating  does  not  depend  on  the  thermal 
conduction  properties  of  the  material,  and  thus,  can  be  particularly  advantageous  in 
processing  low  thermal  conductivity  materials.  The  electromagnetic  power  density 
deposited  in  the  material  is  given  by  Vankoughwett  [8]. 


P  =  2k  f  eQ  e"  E 


Watts 

cubic  meter 


The  effective  dielectric  loss  factor,  e",  a  property  of  the  material,  its 
moisture  content,  and  moisture  distribution,  determines  the  magnitude  of  dielectric 
heating  for  a  given  frequency  and  electromagnetic  field  intensity.  Materials  with  a 
high  affinity  for  electromagnetic  energy  are  called  lossy  materials.  The  e"  depends 
on  the  frequency,  temperature,  moisture  content,  and  electrical  conductivity  of  the 
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Figure  2.2  Schematic  representation  of  the  dependence  of  e"  on  T  and  X 


material,  and  may  also  vary  locally  within  the  material  where  moisture  and 
temperature  gradients  exist  Figure  2.2  shows  a  typical  dependence  of  e"  on  X  and 
T.  For  most  practical  purposes,  an  average  e"  through  the  material  may  be  used 
and  may  be  taken  to  be  a  linear  function  of  X;  consequently,  Qg  may  be 
approximated  as  a  linear  function  of  X,  provided  that  T  and  E  remain  relatively 
constant. 


Heating  of  moist  materials  with  electromagnetic  energy  is  so  effective  that 
reaching  the  boiling  point  is  the  most  common  mode  of  operation  rather  than  the 
exception.  When  the  product  reaches  the  boiling  point,  mass  transfer  of  moisture 
into  the  adjacent  air  is  rapid,  and  the  possibility  of  saturating  the  airstream  becomes 
significant.  Therefore,  the  air  inside  the  dryer  must  be  monitored  if  condensation  is 
to  be  avoided. 

Boiling  point  drying,  characterized  by  high  internal  temperatures  and 
rapid  mass  transfer  due  to  high  energy  input  and  pressure  build-up  in  the  material, 
can  induce  damaging  stresses  in  the  internal  structure  of  the  product  The 
feasibility  of  boiling  point  drying  is  dictated  by  the  product's  thermal  and  physical 
characteristics.  Granular  and  other  porous  materials  are  well-suited  for  boiling 
point  drying  since  internal  pressure  build-up  is  not  severe.  Compacted  materials 
are  more  sensitive  to  internal  pressure  build-up,  thereby  limiting  the  level  of  power 


input.  Materials  that  are  sensitive  to  high  temperatures  may  require  drying  under 
subcooled  conditions  (i.e.,conditions  in  which  the  internal  temperature  is 
intentionally  kept  below  the  boiling  point). 


2.3.3  Dielectrically  Enhanced  Drying  Systems 

A  hybrid  system  is  a  convective  dryer  assisted  by  dielectric  heating.  The 
temperature  difference  between  the  air  and  the  product  will  help  or  hinder  dielectric 
drying.  If  Ta  >  Tp,  evaporation  will  take  place  on  the  surface  and  dielectric  drying 
will  be  enhanced.  The  use  of  dielectric  heating  can  accelerate  the  migration  of  water 
to  the  surface  due  to  internal  pressure  build-up  and  capillary  effects.  Removal  of 
water  by  evaporation  requires  more  energy  than  by  migration  (mechanical  means). 
Taking  into  consideration  the  high  capital  cost  of  dielectric  heating  equipment,  it  is 
generally  less  costly  to  produce  hot  air  than  electromagnetic  energy.  Therefore,  the 
most  efficient  way  to  dry  a  material  is  to  use  dielectric  heating  to  induce  migration 
and  to  use  hot  air  for  evaporation.  Migration  versus  evaporation  is  discussed  by 
Lefeuvre  [10],  but  he  does  not  consider  the  possibility  of  having  a  dry  region 
(evaporative  front)  hindering  the  diffusion  of  heat  and  moisture.  If  a  dry  region 
does  exists,  the  comparison  just  made  between  hot  air  and  dielectric  heating  needs 
to  be  re-evaluated  because  less  heat  from  the  hot  air  is  being  utilized  for  evaporation 
due  to  the  dry  region's  conductive  resistance. 

If  Ta  <  Tp,  condensation  of  the  moisture  vapor  from  the  warm  interior  of 
the  product  will  take  place  at  the  surface  and  a  part  of  the  dielectrically  induced 
heating  will  be  used  to  replace  the  sensible  heat  lost  by  the  material.  There  are 
cases,  however,  in  which  the  cooling  effect  of  the  air  is  desired  to  control  the 
temperature  of  the  product.  The  way  that  air  temperature  is  coupled  with  dielectric 
heating  will  be  an  important  factor  in  the  successful  utilization  of  dielectric  energy 
in  the  drying  industry. 

Schmidt  and  Accad  [9]  have  surveyed  dielectric  drying  modeling  and 
experimental  work.  They  carried  out  parametric  studies  to  explore  the  sensitivity  of 
the  evaporation  rate  from  a  pan  of  water  to  changes  in  Ua,  Ta,  y,  and  Qg.  These 
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Figure  2.3  Relationship  between  Qv,  Qg  and  X. 


studies  showed  that  similar  evaporation  rates  can  be  obtained  for  pure  convective, 
high-temperature  air  drying  and  dielectric  drying  with  ambient  air.  A  “saturation 
point”  was  defined  for  a  dielectric  dryer  when  all  the  incident  power  is  being 
absorbed  by  the  load;  above  the  “saturation”  moisture  content,  Qv  increases  as  the 
water  volume  decreases  while  Qg  stays  approximately  constant  When  the  water 
volume  is  less  than  the  “saturation  point”,  only  part  of  the  incident  power  is 
absorbed  and  the  rest  is  reflected  away  or  absorbed  by  the  surroundings.  At  this 
point,  Qv  stays  constant  (given  that  E  stays  constant),  but  Qg  decreases  as  water 
volume  decreases.  Figure  2.3  presents  a  schematic  representation  of  the 
relationship  between  Qg,  Qv  and  X.  The  value  of  will  later  be  shown  to 
represent  an  important  transition  point  in  the  characterization  of  dielectrically 
enhanced  drying  curves. 

Perkin  [11]  presents  the  continuum  equations  for  mass  flow  of  moisture, 
heat  flow,  and  total  pressure  in  which  the  transport  coefficients  are  a  function  of 
moisture  content  and  temperature  and  must  be  determined  experimentally. 

Perkin  defines  the  parameter  Q^,  which  is  the  equilibrium  temperature 
difference  between  the  product  surface  and  the  air.  Three  drying  regimes  are 
identified  using  this  parameter: 


Qg/h  >  500°C:  Drying  characteristics  are  due  to  electromagnetic  energy  alone,  as 
drying  is  insensitive  to  air  temperature  and  air  humidity.  The  principal  role  of  the 
airstream  is  to  remove  the  evaporating  moisture  from  the  dryer. 


•  Qg/h  <  50°C:  Drying  is  sensitive  to  changes  in  humidity  and  air  temperature. 
Hot  air  contributes  to  the  evaporative  energy. 


•  50°C  <Qg/h  <  500°C:  Behavior  is  in-between  the  above  two  cases. 


Perkin  uses  Bim  to  classify  the  materials  into  either  low  vapor  flow 
internal  resistance  or  high  vapor  flow  internal  -  distance.  The  extreme  cases  for 
Qg/h  and  Bim  are  used  to  simplify  the  governing  equations  of  the  drying  problem 
and  to  obtain  estimates  of  the  drying  behavior  of  the  product 


Pourhiet  and  Bories  [12]  use  a  mathematical  model  originally  established 
by  Luikov  and  experimentally  validated  by  Bories.  The  model  describes  the 
physical  laws  of  heat  and  mass  transfer,  taking  into  consideration  electromagnetic 
power  deposition.  Pourhiet  and  Bories  developed  a  relationship  between  Qg  and  X 
and  speculated  that: 


the  moisture  distribution  inside  the  material  can  be  monitored  using  Qg) 


the  appearance  of  a  dried  zone  near  the  surface  can  be  avoided  by  controlling  Qg, 
and 


•  a  uniform  distribution  of  the  humidity  inside  the  product  can  be  attained  by 
controlling  Qg. 


They  also  found  the  equivalent  airstream  temperature  for  a  given 
microwave  power,  which  would  give  the  same  drying  rate  without  microwaves. 
Therefore,  dielectric  drying  transfer  units  can  be  converted  into  conventional  drying 
transfer  units,  allowing  for  direct  comparison  of  the  two  methods. 


Pourhiet  and  Bories  also  compared  a  hybrid  heating  system  to  a 
conventional  one  to  determine  which  would  be  more  economical  if  both  achieved 
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the  same  drying  time.  Results  showed  that  a  minimum  installed  power  necessary 
to  break  even. 


No  literature  is  currently  available  to  guide  the  development  of 
characteristic  drying  curves  associated  with  dielectrically  assisted  drying.  Statistical 
analysis  techniques,  however,  were  used  by  Meiners  [7]  in  an  attempt  to  normalize 
the  drying  curve  for  a  bed  of  3  mm  glass  beads  under  0.5  kW  of  microwave  power 
and  different  air  temperatures.  While  the  technique  is  useful  in  attempting  to 
extrapolate  data  from  a  limited  series  of  tests  to  a  wider  range  of  conditions,  it  does 
not  provide  the  physical  basis  needed  to  generalize  to  conditions  well  outside  the 
limits  of  available  data. 

The  general  concept  of  the  characteristic  drying  curve  may  apply  to 
dielectric  heating  for  several  reasons: 

•  Under  many  conditions,  the  drying  rate  will  be  dominated  by  the  heat  generation 
term,  which  does  not  depend  on  the  hygrothermal  characteristics  of  air. 

•  Conduction  of  convective  heat  into  the  material  will,  in  many  cases,  not  be  the 
primary  driving  force  for  evaporation;  therefore,  the  depth  of  the  dry  region, 
which  is  established  by  Q,  will  not  be  the  thermally  limiting  element 

Keey  [6,16]  considered  these  parameters  as  impediments  to  the 
development  of  a  single  characteristic  drying  curve  for  a  particular  product  In 
dielectrically  enhanced  drying,  the  dielectric  response  characteristic  of  the  material 
and  the  interrelationship  between  internal  and  surface  heating  (or  cooling)  are 
expected  to  play  the  primary  role  in  defining  drying  behavior.  These  parameters 
will  be  considered  in  detail  in  Chapter  3. 


CHAPTER  3 


CHARACTERISTIC  DRYING  CURVE  FOR  DIELECTRICALLY 

ENHANCED  DRYING 

3.1  NONDIMENSIONALIZATION  OF  THE  PROBLEM. 

Proper  selection  of  normalizing  scales  produces  variables  of  comparable 
magnitude  over  a  range  of  operating  conditions.  The  pertinent  parameters  for 
nondimensionalization  of  a  problem  can  be  found  either  by  bringing  into  play 
information  provided  by  basic  physical  laws,  or  assumptions  about  the  functional 
form  of  the  result  [18].  The  objectives  of  nondimensionalizing  the  problem  are  to 
reduce  the  number  of  variables,  to  simplify  the  mathematical  structure  of  the 
problem  to  make  results  more  widely  applicable,  and  to  keep  experimental  runs  to  a 
minimum.  In  drying  processes,  the  main  motivation  for  this  approach  is  to  limit  the 
amount  of  experimental  data  needed  to  extrapolate  results  from  a  particular  drying 
condition  to  other  drying  condition'.  Eckert  [19]  discusses  the 
nondimensionalization  of  a  convective  drying  problem  in  which  the  difference 
between  the  wet  bulb  temperature  and  the  product's  initial  temperature  is  used  to 
make  the  conservation  equations  dimensionless.  The  nondimensionalization  of  the 
problem  allowed  Eckert  to  perform  parametric  studies  by  varying  a  few 
dimensionless  parameters  instead  of  the  numerous  variables  in  the  problem. 

In  this  chapter,  an  energy  balance,  performed  on  a  control  volume  element 
representative  of  the  wet  product  with  convection  and  dielectric  heating,  is  used  to 
develop  dimensionless  parameters  with  the  intention  of  using  these  parameters  to 
produce  a  characteristic  drying  curve. 

Refering  to  Figure  3.1,  an  energy  balance  on  the  control  volume  around 
the  wet  product  yields 
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where  all  the  fluxes  are  based  on  the  evaporative  surface  area,  and 
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Qw  =  psb(Cs  +  XCw)-^  (3.4), (3.5) 


In  Equation  3.1,  the  change  of  energy  content  associated  with  water 
vapor,  air,  and  phase  transformation  of  the  solid  species  has  been  neglected-  Note 
that  in  Equation  3.3,  Qg  was  defined  in  terms  of  the  dielectric  forward  power  which 
is  assumed  to  be  constant  and  completely  absorbed  by  the  moist  product  in  the 
initial  drying  stage;  for  later  stages  where  Qg  is  expected  to  decrease,  the 
characteristic  drying  curve  concept  will  be  used  to  represent  the  decline  in  Qg 
departing  from  the  maximum  value  calculated  in  Equation  3.3.  Some  of  the 
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variables  in  Equation  3.1  are  hard  to  normalize  because  they  do  not  have  an 
obvious  scaling  factor.  An  evident  choice  for  a  normalizing  parameter  would  be  a 
direct  operating  condition  since  it  is  known  and  can  usually  be  set  to  a  constant 
value.  Examples  of  such  operating  conditions  are  air  temperature,  air  humidity, 
microwave  or  radio  frequency  foward  power,  air  velocity,  and  initial  moisture 
content 

The  characteristic  drying  curve  is  defined  in  terms  of  Nv  and  X,  as 
discussed  in  Chapter  2.  For  convenience,  therefore.  Equations  3. 1-3.5  will  be 
transformed  from  their  independent  variable,  t,  to  the  more  significant  variable,  X. 
To  transform  from  “time  domain”  to  “moisture  content  domain”: 


Therefore,  Equation  3.5  can  now  be  expressed  in  terms  of  X; 


Q  =  -N  (C  +  X  C  ) 
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Combining  Equations  3.7,  3.4  and  3.1, 
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The  LHS  of  Equation  3.8  is  the  total  energy,  Qr,  going  Into  the  product.  The  RHS 
of  Equation  3.8  will  depend  on  Qr  regardless  of  the  heat  source,  dielectric  or 
convective.  Qr  could  be  dominated  by  Qc  or  Qg;  in  this  case,  Qr  will  be 
normalized  on  the  dielectric  heat  flux. 


Note  that  if  Qt  is  dielectrically  dominated,  0-1,  and  if  Qt  is  convectively 
dominated,  0  »  2.  This  normalized  parameter  will  be  termed  the  “drying  mode 
factor”.  Its  physical  significance  in  dielectrically  assisted  drying  will  be  discussed 
at  length  in  Section  3.2. 

Using  Equation  3.9  to  replace  (Qc  +  Qg)  and  dividing  Equation  3.8  by 

0Qg: 
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(3.10 


The  quantity  within  the  squared  brackets  is  the  specific  energy  (kJ/kg)  gained  by 
the  wet  product  normalized  by  the  latent  heat  of  vaporization  of  water.  This  term 
includes  a  temperature  gradient  characteristic  of  the  warm-up  period.  Define  the 
bracketed  term  as. 
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The  dimensionless  variable,  fw,  will  be  called  the  “relative  warm-up  rate”.  This 
parameter  is  not  found  in  the  drying  literature;  its  physical  significance  in  dielectric 
drying  is  discussed  in  detail  in  Section  3.1.1. 


Going  back  to  Equation  3.10,  the  first  term  in  the  RHS  of  the  equation  is 
the  drying  rate,  Nv,  normalized  by  the  maximum  drying  rate  when  the  moisture  in 
the  solid  reaches  saturation  temperature.  Define  the  dimensionless  variable,  /r,  as: 
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/r ,  will  be  called  the  “relative  drying  rate”,  which  behaves  similarly  to  its 
counterpart  in  convective  drying.  Finally,  the  energy  balance  equation,  3.1,  takes 
the  normalized  form: 


1  =  V'-/w> 


(3.13) 


The  following  sections  explore  the  relationship  of  the  normalized  parameters,  /w, 
and  /r,  to  the  characteristic  drying  curve  concept. 


To  study  the  physical  meaning  of  /w  ,  the  variables  Tp  and  X  will  be 
nondimensionalized  using  the  range  of  these  variables  during  the  drying  process  as 
scaling  factors.  For  Tp,  the  scaling  variable  is  the  water's  local  saturation 
temperature,  Tb,  and  the  reference  variable  is  the  initial  product  temperature,  T0. 
For  X,  scaling  is  done  by  the  boiling  point  moisture  content,  Xb,  and  the  reference 
value  is  given  by  the  initial  moisture  content,  Xq.  That  is. 


*  T  -T 
X  = _ E _ JL 

P  TD  -T 
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*  X-X 
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Performing  a  variable  transformation  from  the  pair  (X,  T)  to  the  pair  (X*,T*)  the 
following  expression  is  obtained  for  /w: 


Av  = 
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(3.15) 


Two  groups  of  variables  generated  by  the  normalization  are  further  defined: 
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Now  Equation  3.15  reads: 
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(3. 16),(3. 17) 


/w  =  Q*  (  C*  +  X* )  — E-  (3.18) 

ax 

The  quantity  Q*  is  the  dimensionless  ratio  of  the  average  sensible  heat  gained  by  the 
water  inside  the  material  from  t  =  0  to  saturation  to  the  latent  heat  of  vaporization. 

C*  is  the  ratio  of  the  initial  thermal  mass  of  the  wet  product  to  the  total  thermal  mass 
of  the  evaporated  water  from  t  =  0  to  saturation. 

The  normalized  temperature  gradient  in  Equation  3.18  will  now  be 
examined.  From  a  purely  intuitive  point  of  view,  it  can  be  argued  that  if  Qr 
increases,  ATp  will  increase,  but  AX  will  also  increase,  making  their  ratio  constant 
for  a  given  material  at  a  given  moisture  content  Therefore,  3T/9X  is  only  a 
function  of  X  and  the  properties  of  the  material.  Furthermore,  Q*  and  C*  depend 
only  on  the  material  properties  and  the  initial  conditions .  From  these  observations, 
the  following  assertion  is  made:  for  a  given  material,  the  relative  warm-up  rate  is 
independent  of  air  temperature  and  dielectric  power  and  only  depends  on  the 
moisture  content.  To  test  this  hypothesis,  9T* /3X*  vs.  X  and  /w  vs.  X  were 
plotted  in  Figures  3.2  and  3.3  respectively,  using  drying  data  for  3  mm  glass  beads 
from  several  experiments  run  at  different  air  temperatures  and  power  densities,  with 
T0  »  25°C.  As  expected,  dT* /9X*  and  /w  behaved  in  the  same  fashion  regardless 
of  the  power  density  or  air  temperature.  Therefore,  at  least  for  the  limited  set  of 
data  given  here,  the  relative  warm-up  rate  is  a  function  only  of  X  during  the  warm¬ 
up  period. 


ll 


The  normalized  drying  rate,  /r,  called  the  relative  drying  rate,  approaches 
one  as  the  relative  warm-up  rate  approaches  zero.  At  this  point,  in  boiling  point 
drying,  the  drying  rate  is  observed  to  become  constant.  Figure  3.4  was  generated 
using  drying  data  for  a  bed  of  glass  beads  0.1  m2  in  surface  area  and  5  cm  deep, 
subject  to  0.5  kW  microwave  power  and  Ta  =  93°C.  Note  that  the  constant  drying 
rate  begins  approximately  at  the  point  where  the  wet  product  reaches  saturation 


temperature. 


At  a  certain  moisture  content,  all  the  forward  dielectric  power  cannot  be 
absorbed  by  the  sample,  and  some  of  the  power  is  reflected.  The  generated  power 
starts  decreasing  as  e"  decreases  with  decreasing  X  (see  discussion  in  Chapter  2). 
At  this  point,  Nv  starts  decreasing  proportionally  to  Qg  since  no  sensible  heating  is 
occurring.  This  period  of  the  dielectric  drying  process  is  referred  to  as  the  falling 
rate  period.  The  falling  rate  period  in  this  case  is  produced  primarily  by  a  decline  in 
Qg  and  not  by  limitations  on  moisture  diffusion  in  the  product  The  moisture 
content  at  which  the  power  deposition  starts  decreasing  is  here  called  the  “critical 
moisture  content”,  Xcd.  It  has  the  same  characteristics  as  the  critical  moisture 
content  in  conventional  drying,  but  the  mechanisms  behind  them  are  completely 
different  It  can  be  inferred  that  during  the  falling  rate  period,  /r  is  only  a  function 
of  X. 


In  summary,  Sections  3.1.1  and  3.1.2  have  asserted  that  fT  is  a  function 
only  of  X  throughout  the  drying  process  which  implies  that  there  is  a  unique 
characteristic  drying  curve  for  materials  dried  under  dielectric  power.  Experimental 
evidence  was  presented  for  a  typical  nonhygroscopic  material.  Going  back  to 
Equation  3. 13,  /r  will  be  a  maximum  when  /w  =  0;  then,  fT  =  1  until  X  =  Xcd-  At 
this  point,  the  power  input  is  no  longer  constant  and  /r  decreases  until  X  =  0. 

Since  it  was  demonstrated  that  /r  and  /w  only  depend  on  X  for  a  given  material,  the 
following  expression  can  be  stated: 
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(3.19) 


where  91  (X)  is  obtained  by  means  of  curve  fitting  or  estimating  techniques 
performed  on  experimental  test  data  as  discussed  in  Section  3.3.  Note  that  91  (X) 
goes  from  =  0  to  1  in  the  warm-up  period;  is  equal  to  one  in  the  constant  rate 
period;  and  goes  from  1  to  ~  0  in  the  falling  rate  period. 


3.2.1  Physical  significance 

The  drying  mode  factor  indicates  whether  drying  is  dominated  by 
convective  or  dielectric  heating.  This  parameter  is  also  used  as  a  scale  to  normalize 
different  convective  and  conductive  conditions,  and  to  define  different  drying 
regimes  that  arise  as  a  consequence  of  dielectric  heating.  Also,  0  is  used  to  estimate 
the  temperature  profile  in  the  product's  dry  region.  This  particular  effect  is 
explored  in  Section  4.1.3. 

Note  that  0  is  the  normalized  Qr  going  into  the  wet  product  which  is 
influenced  by  the  presence  or  absence  of  an  evaporative  front  due  to  the  conductive 
resistance  imposed  by  the  dry  region.  If  a  dry  region  is  not  present  then  Equation 
3.9  becomes 


0=1  + 


h(Ta-T) 


(3.19) 


However,  if  a  dry  region  develops  inside  the  material,  the  expression  for  Qr  will 
change  because  the  convective  heat  flux  is  now  being  used  for  both  sensible  heating 
in  the  dry  region  and  sensible  heating  and  evaporation  in  the  wet  region. 

Therefore,  if  a  dry  region  is  present. 


Figure  3.5  Control  volume  for  the 


Note  that  <(>  =  0'  when  an  evaporative  front  is  not  present  since  Qc  =  Q§.  Also,  <j) 
and  <J>'  approach  infinity  as  Qg  approaches  zero.  This  singularity  does  not  present  a 
problem  because  as  <}>  approaches  infinity,  the  product  <}>  Qg  approaches  Qc  in 
Equation  3.10. 

To  calculate  <J>’,  it  is  necessary  to  know  the  instantaneous  position  of  the 
evaporative  front  and  the  temperature  profile  inside  the  dry  region.  If  this 
information  is  not  available,  and  the  thermal  mass  of  the  dry  product  is  relatively 
small,  the  convective  heat  flux  can  be  used  as  an  approximation  for  the  conducted 
heat  without  incurring  a  large  error.  The  energy  balance  on  the  dry  region  will  be 
used  to  establish  the  relationship  between  <f>  and  <j>'.  From  Figure  3.5: 


Equation  3.22  gives  the  sensible  heat  gained  by  the  dry  region  based  on  the  mean 
temperature  of  the  region;  the  thermal  effect  of  the  region's  growth  was  not 
included.  Transforming  from  time  domain  to  moisture  content  domain. 


3T  -N  3T 
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Combining  Equations  3.21,  3.22,  and  3.23,  dividing  by  Qg,  adding  one  to  both 
sides  and  using  Equation  3.19  and  3.20,  the  following  expression  is  obtained: 
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From  equation  3.24,  <j)  =  <f>'  as  long  as: 

•  aTm/ax  ~  0,  which  translates  into  negligible  sensible  heating  in  the  dry  region 

•  5  /  b  is  very  small,  or 


•  Qg  is  much  bigger  than  Qc  (<$)  =  1) 

A  drying  mode  factor  cannot  be  calculated  for  cases  in  which  the  dry 
region  thickness  is  large  and  unknown.  For  such  cases,  a  reasonable  assumption 
would  be  <b'  =»  1  since  a  large  conductive  resistance  is  expected,  making  the 
conduction  heat  transfer  small.  Nou.  that  the  drying  mode  factor  is  not  necessarily 


constant,  since  it  depends  on  Ts.  Product  temperatures  need  to  be  measured  and  an 
empirical  or  analytical  expression  is  needed  to  estimate  the  position  of  the 
evaporative  front. 


3.2.2  Drying  Regimes 

Perkin  [11]  discusses  the  usefulness  of  defining  drying  regimes  for 
dielectrically  assisted  drying.  To  define  the  regimes,  he  uses  the  ratio  Qg/h.  To 
date,  Perkin  s  approach  has  been  the  only  effort  to  quantify  the  importance  of  the 
complimentary  relationship  between  dielectric  heating  and  convective  heating  in 
terms  of  efficient  drying.  The  disadvantage  of  the  Qg/h  approach  is  that  it  is  not 
dimensionless  and  does  not  take  into  account  the  relationship  between  the  product 
and  air  temperatures. 

The  drying  mode  factor  will  be  used  to  define  distinct  drying  regimes  that 
arise  from  dielectric  drying.  The  following  is  a  compilation  of  limiting  cases  and 
the  significance  of  the  drying  mode  factor  for  each  regime. 

•  Convective  Regime  (0  »2):  The  product  experiences  a  high  rate  of  convective 
heating  and  the  convective  heat  flux  is  much  higher  than  the  dielectric  heat  flux. 
This  regime  corresponds  to  Perkin's  drying  regime  defined  by  Qg/h  <  50°C. 
Drying  is  sensitive  to  changes  in  air  conditions. 

•  Mixed  Regime  (0  =  21:  Convective  and  dielectric  heating  are  approximately 
equal  in  magnitude. 

•  Dielectric  Regime  (0=1):  This  regime  is  characterized  by  high  dielectric 
heating  (Qg  »  Qc)  or  negligible  convective  heating.  This  corresponds  to  Perkin’s 
drying  regime  defined  by  Qg/h  >  500°C.  Drying  is  relatively  insensitive  to  air 
conditions,  and  the  principal  function  of  the  air  is  to  serve  as  a  moisture  sink. 

•  Subcooled  Regime  (  0  <  1  ):  The  product  experiences  convective  cooling.  If 
0  =  0,  the  convective  cooling  and  dielectric  heating  are  approximately  equal  in 
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magnitude.  The  product  temperature  is  higher  than  the  dry  bulb  temperature  of 
the  air  which  can  only  occur  when  dielectric  heating  is  applied 

•  Cooling  Regime  (<j>  <  0):  This  implies  that  the  product  temperature  is  decreasing 
as  it  is  being  dried.  An  example  of  this  regime  is  the  cool-down  zone  in  a  multi¬ 
zone  dryer,  where  the  product  is  allowed  to  approach  ambient  temperature, 
eliminating  remaining  moisture  without  dielectric  heating. 

H  PROCEDURE  FOR  NORMALIZING  EXPERIMENTAL  DATA 

3.3.1  Collection  of  experimental  data 

The  experimental  data  required  are  mw,  Ts,  Ta,  Ua,  and  P,  as  a  function 
of  t.  The  data  are  recorded  using  a  computerized  data  acquisition  system.  After 
placing  the  wet  sample  on  top  of  the  digital  balance,  the  balance  is  tared  to  monitor 
the  amount  of  water  that  has  evaporated.  Temperature  probes,  connected  to  a 
Luxtron  fiber-optic  temperature  measurement  system,  are  placed  on  the  surface  and 
in  the  middle  of  the  wet  sample.  A  velocity  probe  is  used  to  monitor  the  air 
velocity,  and  a  dual  channel  microwave  power  meter  is  used  to  monitor  the  forward 
and  reflected  power  going  into  the  dryer.  The  data  are  written  to  magnetic  media  at 
predetermined  time  intervals. 

3.3.2  Data  transformation. 

The  experimental  data  must  be  transformed  into  X  and  Nv  terms.  To 

obtain  X, 


X  = 
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where  mj  is  the  mass  of  the  evaporated  water  at  tj.  To  obtain  Nv, 
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A  smoothing  technique  is  used  to  reduce  extrinsic  noise  inherent  in  the  experimental 
procedure.  An  average  value  from  n  consecutive  data  points  is  used  to  represent  an 
interval. 
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For  example  if  n  =  5  and  i  =  7  then 
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Finally,  a  data  file  is  created  containing  X,  Nv,  Ts,  and  Tp,  for  a  particular  Ta  and 
P. 

3.3.3  Data  normalization 

The  data  file  is  read  by  the  normalizing  code,  and  the  normalizing 
parameters  are  calculated  for  each  data  point.  Letting  Qg  =  P/A, 
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where  hj  is  calculated  by  a  subroutine  based  on  empirical  correlations  of  the  Nusselt 
number  given  in  the  form 


Nu  =  Cj  Re  Pr 


Finally,  the  empirical  data  are  normalized  by  the  expression, 
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The  code  used  to  normalize  the  data  is  listed  in  Appendix  A. 


The  normalized  data  must  be  plotted  to  assess  the  accuracy  of  the 
normalization  and  to  obtain  the  characteristic  drying  curve  equation.  Few  the  warm 
up  period,  a  characteristic  function  is  determined  from  Xo  to  Xb-  Typically  a 
linear,  parabolic  or  cubic  function  is  obtained;  so,  fT  can  be  represented  as: 


=  (1'4)X  *f0 
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fr  =  (l-/o)X  -3(l-/o)X  +3(l-/o)X  +/o  (3.34 


where  f0  is  the  initial  relative  drying  rate.  Functions  3.32  through  3.34  are 
graphed  in  Figure  3.6,  where  it  was  assumed  that  f0  =  0.  Note  that  X*  is  only 
defined  for  the  warm-up  period. 

For  the  falling  rate  period,  Accad  and  Schmidt  [9]  and  others  (Figure  3.4) 
have  observed  that  the  drying  rate  decreases  linearly  with  decreasing  X  (See 
Chapter  2);  however,  cubic  functions  have  also  been  obtained  for  the  falling  rate 
period  (see  Section  5.3.3).  The  characteristic  drying  curve  function  can  be 
obtained  by  applying  either  equation  3.32,  3.33,  or  3.34  and  by  substituting  f0  by 
/f.  The  normalized  X  for  the  falling  rate  period  is  given  by 

X  -X 

O  =  — i-  (3.35) 
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Note  that  for  Xf  =  0,  O  becomes  the  normalized  moisture  content  as  defined  for 
convective  drying  of  nonhygroscopic  materials. 
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To  illustrate  the  proposed  normalization  scheme,  existing  sets  of 
experimental  data  on  a  bed  of  3  mm  glass  beads  0.1  m2  in  area  and  5  cm  deep, 
subject  to  different  microwave  power  levels  and  air  temperatures,  are  normalized. 
The  experimental  data  are  limited  in  the  sense  that  it  only  covers  the  range  of 
0.8  <  <t>  <  1.25,  but  it  does  demonstrate  the  normalization  scheme  for  different 
power  levels  and  air  conditions.  Since  the  bed  of  glass  beads  is  a  coarse 
nonhygroscopic  material,  an  evaporative  front  is  expected  to  develop  inside  the 
product  The  function  proposed  in  Equation  4.10  is  used  to  estimate  the  front 
position,  with  the  evaporative  front  coefficient,  (3,  estimated  on  the  basis  of  air 
temperature  and  delivered  microwave  power  (see  Section  4.1.1).  Evidence  of  the 
existence  of  an  evaporative  front  in  the  material  can  be  observed  in  Figure  3.7 
where  the  drying  rate  curves  for  different  temperatures  and  power  densities  have 
been  plotted.  Note  that  for  P  =  0.5  kW,  the  drying  curves  start  together,  grow 
apart  from  each  other  in  proportion  to  (Ta  -  Ts),  and  finally  converge  into  one 


again.  If  a  drying  test  could  be  set  up  where  the  air  temperature  is  always  equal  to 
the  product  temperature,  the  convective  heat  flux  would  be  zero  at  all  times.  Under 
these  conditions,  the  drying  rate  curve  is  expected  to  rise  until  no  sensible  heating  is 
occurring.  At  this  point,  all  the  energy  input  goes  into  evaporation,  and  the  drying 
rate  will  be  constant.  However,  heating  or  cooling  effects  will  cause  the  drying  rate 
to  depart  from  the  pure  dielectric  drying  behavior.  If  convective  cooling  is  present, 
the  energy  left  for  evaporation  is  less  than  the  heat  generated  and  the  drying  rate 
will  be  proportionally  lower,  if  convective  heating  is  present,  the  energy  available 
for  evaporation  will  be  more  than  the  heat  generated,  and  the  drying  rate  will  be 
proportionally  higher.  This  observation  can  be  seen  in  Figure  3.7.  The  curves 
grow  in  proportion  to  the  heat  input  until  a  moisture  content  of  about  15%  is 
reached.  After  this  point,  the  curves  tend  to  converge  and  follow  the  same  path 
regardless  of  the  air  temperature;  therefore,  Qc  is  not  contributing  to  Qr-  At  the 
beginning  of  the  drying  period,  the  dry  region  is  relatively  small  and  the  heat 
conduction  into  the  material  is  not  hindered-  At  some  point,  the  dry  region 
becomes  large  enough  to  impede  the  conduction  of  heat  to  or  from  the  wet  region, 
and  dielectric  heating  is  practically  the  only  source  of  evaporation.  This  situation 
resembles  the  hypothetical  case  of  pure  dielectric  drying  mentioned  earlier  and  all 
the  drying  curves  will  converge  into  one. 

The  data  collected  for  the  drying  tests  were  transformed  from  the 
parameters  (t,mw,T)  to  the  parameters  (X,NV,T)  and  input  into  the  normalizing 
code.  Figure  3.8  shows  the  results  of  the  normalization.  Despite  the  experimental 
uncertainties  of  the  measuring  devices,  the  warm-up  rate  period  normalized  fairly 
well.  A  constant  drying  rate  period  was  not  observed  for  the  0.5  kW  curves 
whereas  a  well-defined  constant  drying  rate  was  observed  for  the  1.0  kW  curves. 

The  results  shown  in  Figures  3.7  and  3.8  support  the  hypothesis  that  the 
normalization  procedure  outlined  in  this  chapter  can  be  used  to  characterize  drying 
rates  for  combined  dielectric  heating  and  convective  conditions. 
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Figure  4. 1  Illustration  of  the  receding  plane  model. 


The  following  equation  relates  the  drying  rate  to  the  product’s  physical 
dimensions  and  the  main  independent  variable,  X: 


0> 


Equation  4.2  relates  the  drying  rate  to  the  physical  dimensions  of  the  wet  product. 


Note  that  pw  is  defined  as  the  mass  of  water  per  unit  volume  of  the  wet  region  (not 
total  product  volume),  and  hence,  depends  on  how  the  moisture  is  distributed  in  the 
product  Consider  the  two  extreme  cases: 

Case  1.  Constant  water  density.  This  case  is  characterized  by  a  well-defined 
evaporative  front  sweeping  down  the  product  with  the  wet  region  pore  volume 
totally  saturated.  Equation  4.2  reduces  to: 


Combining  Equation  4. 1  and  4.3  and  solving  for  85: 
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The  ratio  of  the  densities  in  Equation  4.4  is  given  by: 
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Now,  Equation  4.4  can  be  arranged  into  an  integrable  form. 
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Integrating  Equation  4.6  yields: 


£-=  5*  =  ( 1  -  ) 
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Equation  4.7  states  that  the  evaporative  front  position  is  a  linear  function  of  the 
moisture  content  for  this  particular  case. 

Case  2.  Volume  of  the  wet  region  remains  constant  The  moisture  content  is 
spatially  uniform;  in  other  words, 


ax  =  ax  _ax 

8x  8y  8z 


(4.8) 
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This  particular  case  was  observed  by  Lyons  [13]  when  he  subjected  a  bundle  of 
wet  cotton  yam  to  microwave  heating.  Throughout  the  experiment,  the  moisture 
distribution  was  uniform  inside  the  material  (See  Figure  4.2).  For  this  case. 
Equation  4.2  reduces  to 


dp 

AN  =  -V 
v  w  ^ 


(4.9) 


and  going  through  the  same  procedure,  8*=  0  (no  evaporative  front  is  present). 

Most  materials  are  expected  to  behave  somewhere  between  Case  1  and 
Case  2  depending  on  their  permeability,  dielectric  heating  rate  and  air  temperature. 
Low  permeability  and  high  power  will  generate  pressure  inside  the  material  which 
will  tend  to  keep  the  surface  wet  and  8  small.  Also,  if  boiling  is  occurring  within 
the  wet  region,  the  boiling  action  tends  to  push  the  water  towards  the  surface.  This 
behavior  has  been  observed  experimentally  in  drying  of  beds  of  glass  beads  under 
microwave  irradiation  [20].  Also  for  low  surface  temperature,  continuous  vapor 
condensation  in  the  dry  region  will  keep  this  region  wet,  raising  the  front  posidon. 
In  addition,  if  the  permeability  of  the  material  is  low,  pressure  build-up  and 
capillary  forces  will  drive  moisture  toward  the  surface,  hindering  the  establishment 
of  an  evaporative  front. 

Based  on  the  above  observations,  the  following  model  is  proposed: 

cj*  J  /  X  \P 

5  =  1  '  x~'  (4.10) 

o 


where  P,  the  evaporative  front  coefficient,  is  a  function  of  the  permeability, 
surface  temperature  and  microwave  power.  In  general: 


P  — » 1  for  granular  nonhygroscopic  materials  with  a  high  surface  temperature; 
P  — >  0  for  low  permeability  materials  and  low  surface  temperature. 
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Figure  4.3  shows  the  position  of  the  evaporative  front  measured  from  the  surface. 
If  5*  =  1,  the  entire  product  is  dry;  if  8*  =  0,  the  moist  region  encompasses  the 
whole  product  volume. 


For  convective  drying,  an  equation  for  8*  is  given  by  Suzuki  [  14]  as: 


8=1- 


x/t  “■[#■] 


(1  -  n) 
2 


The  relationship  fT  =  d>n  was  used  to  develop  Equation  4.11  in  terms  of  X, 
assuming  constant  moisture  diffusivity  in  a  porous  nonhygroscopic  material.  Note 
that  Equations  4.10  and  4.1 1,  are  similar  in  nature.  Evaluation  of  p  for  different 
types  of  materials  as  a  function  of  their  moisture  transport  characteristics  is  not 
within  the  scope  of  this  report.  Efforts  are  presently  being  devoted  by  others  to 
developing  models  for  drying  kinetics  of  solids  which  will  permit  direct  estimation 
of  P  from  fundamental  material  characteristics  [23]. 


The  dielectrically  assisted,  convective  tunnel  dryer  was  divided  into  three 
control  volumes  as  depicted  in  Figure  4.4.  Control  Volume  I  encompasses  the 
humid  air  inside  the  dryer,  Control  Volume  II  corresponds  to  the  dry  region  left 
behind  by  an  evaporative  front  moving  into  the  product,  and  Control  Volume  ID  is 
the  wet  region  of  the  product.  A  lumped  parameter  analysis  was  made  for  the  y 
direction  thermal  and  mass  transport  between  the  regions  [2],  Due  to  the  small 
temperature  gradients,  heat  conduction  in  the  z  and  x  directions  was  neglected  in  the 
analysis.  Equations  4. 1 2,  4. 1 3,  and  4. 1 4  are  the  energy  balance  equations  for 
Control  Volumes  I,  II,  and  HI,  respectively.  The  complete  derivations  of  these 
equations  is  found  in  Appendix  B. 
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where  £  =  1  -  5*.  Equation  4. 13  incorporates  two  new  terms  into  the  existing  dryer 
model.  NVCV  (Tm  -  Tp )  is  the  change  in  enthalpy  experienced  by  the  water  vapor 
while  it  is  being  transported  from  the  evaporative  interface  to  the  surface.  This  term 
was  initially  neglected  by  Flake  [2]  because  Nv  and  (Tm  -  Tp )  are  relatively  small 
in  convective  drying.  The  incorporation  of  dielectric  heat  generates  a  larger  Nv, 
and  the  temperature  difference  (Tm  -  Tp )  is  larger  when  Ta  »  Tb  or  Ta  «  Tg.  These 
two  factors  make  the  change  in  the  vapor  enthalpy  term  an  influential  one. 

The  third  term  in  the  RHS  of  Equation  4. 13  accounts  for  the  dry  product 
enthalpy  influx  into  the  control  volume  due  to  the  moving  evaporative  plane.  This 
term  was  also  neglected  by  Flake  [2]  because  it  is  proportional  to  Nv  and  Tp  which 
are  relatively  small  for  convective  drying.  Examining  this  term  for  the  case  where 
3  =  1,  the  differential  term  is 
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Therefore, 
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Equation  4.16  reveals  that  the  third  term  in  Equation  4.1 1  is  most  significant  for 
less  porous  and  denser  materials  when  experiencing  relatively  high  drying  rates  at 
relatively  high  temperatures. 

Finally,  Equation  4. 14  states  that  the  heat  conducted  into  the  wet  region 
plus  the  heat  generated  must  equal  the  evaporative  heat  loss  plus  the  energy 
associated  with  the  dry  product  leaving  the  control  volume  due  to  the  movement  of 
the  evaporative  front  plus  the  change  in  sensible  heat  of  the  wet  region. 

From  a  composite  mass  balance  on  Control  Volumes  I,  II,  and  III 


V  a  dz  s  dz 


(4.17) 


The  energy  and  mass  balance  relations  expressed  in  Equations  4.12  through  4.14 
and  Equation  4.17  must  be  supplemented  by  a  number  of  additional  constitutive 
relations  to  complete  the  formulation  of  the  problem.  The  equations  to  complete  the 
analysis  were  provided  by  Flake  [2]  and  are  summarized  in  Appendix  B. 


4.1.3  Temperature  profile 


Although  drying  is  a  transient  phenomenon,  quasi-steady  analyses  may  be 
valid  if  the  time  scale  for  the  rate  of  recession  is  large  compared  to  the  time  required 
for  a  steady  state  response  of  the  temperature  and  humidity  gradients  to  the  moving 
boundary  [2].  This  is  generally  valid  for  convective  drying  in  which  the  drying 
time  is  relatively  long,  the  air  temperature  is  higher  than  the  product  temperature, 
and  the  effect  of  vapor  convection  in  the  dry  region's  temperature  profile  can  be 
reasonably  neglected.  When  dealing  with  dielectrically  enhanced  drying,  there  will 
be  many  cases  where  the  product  temperature  is  higher  than  the  air  temperature 
(subcooled  drying),  vapor  fluxes  are  much  larger  than  in  convective  drying,  and 


drying  rate  is  much  higher.  It  has  been  observed  in  recent  experiments  that  the 
vapor  flow  through  the  dry  region  tends  to  shift  the  temperature  profile  towards  a 
more  uniform  state,  depending  on  the  air  temperature  and  dielectric  power.  Figure 
4.5  shows  the  temperature  history  for  drying  of  glass  beads  at  Ta  =  65°C  for 
1.0  kW.  The  equilibrium  temperature  for  the  surface  was  around  80°C  even 
though  Ta  =  65°C.  The  conduction  heat  flux  from  the  evaporative  plane  to  the 
surface  is  not  expected  to  have  a  large  effect  on  the  surface  temperature  because  of 
the  very  low  thermal  conductivity  of  the  glass  beads  (=  0.0004  kW/m-K). 
Therefore,  it  is  concluded  that  for  these  two  cases,  Ts  >  Ta  and  Ts  approaches  Tp 
because  the  vapor  flux  through  the  dry  region  keeps  this  region  at  a  more  uniform 
temperature.  This  behavior  depends  on  Ta  and  <j>;  their  effect  on  the  assumed  linear 
temperature  profile  inside  the  dry  region  can  be  stated  as  follows: 

•  If  Ta  »  Tb  and  <J>  »  2,  the  temperature  profile  is  not  significantly  affected,  and  the 
dry  region  temperature  profile  may  be  assumed  linear. 

•  If  Ta  »  Tb  and  <j)  *  1 ,  the  vapor  flux  tends  to  cool  down  the  dry  region  of  the 
product  and  keep  it  near  saturation  temperature.  In  this  case,  it  is  reasonable  to 
assume  that  the  wet  and  dry  region  are  approximately  at  the  same  temperature. 

•  If  Ta  =  Tb,  <f>  becomes  unity  automatically,  and  the  temperature  profile  throughout 
the  product  is  uniform. 

•  If  Ta  «  Tb,  and  0  «  1,  the  temperature  profile  is  not  affected,  and  the  temperature 
inside  the  dry  region  is  not  uniform. 

•  If  Ta  «  Tb  and  <J>  =  1,  the  vapor  transport  through  the  region  will  tend  to  heat  up 
the  dry  region  and  establish  a  uniform  temperature  throughout  the  product 


In  summary,  when  <j>  *  1,  it  can  be  assumed  that  the  dry  region  temperature  equals 
that  of  the  wet  region. 

Flake  [2]  assumes  that  the  dry  region  is  at  a  mean  temperature  Tm, 
derived  from  the  arithmetic  average  of  Ts  and  Tp.  This  assumption  is  suitable  for 
low  intensity  convective  drying  where  a  front  is  observed  only  in  the  late  stages  of 
drying,  and  any  error  due  to  the  assumption  of  a  linear  temperature  profile  in  the 
dry  region  is  negligible.  However,  for  coarse  nonhygroscopic  materials  under 
dielectric  heating,  a  front  may  occur  from  the  beginning  and  the  dry  product 
entering  the  control  volume  in  the  y  direction  can  have  a  significantly  higher 
enthalpy  than  the  rest  of  the  dry  region.  This  observation  has  been  incorporated 
into  the  governing  equations  of  the  product  model.  Appendix  B  contains  a  more 
detailed  derivation  of  the  analytical  expression  to  account  for  this  "extra”  enthalpy 
carried  into  the  control  volume  by  the  dry  product. 

4.2  HEAT  GENERATION  TERM 

Schmidt  and  Acaad  [9]  performed  a  series  of  model  studies  which 
considered  a  pan  of  water  subjected  to  simultaneous  dielectric  heating  and 
convection.  The  volumetric  heat  generation  term  was  plotted  against  the  water 
volume,  and  a  quasi-linear  dependence  between  the  heat  generation  and  the  water 
volume  was  observed,  as  discussed  in  Chapter  2.  From  this  behavior,  it  is 
assumed,  that  if  X  =  f(z)  then  Qv  =  Tf(z),  where 


Qv  _  P  _  1  3P 
X  =  = 


(4.18) 


The  initial  conditions  are  used  to  determine  T.  Since  T  is  constant,  it  also 
represents  the  ratio  of  Q V/X  for  a  differential  control  volume  element  as  expressed 
by  the  last  term  of  Equation  4.18.  Solving  this  equation  fordP: 


It  is  desirable  to  express  the  rate  of  heat  generated  on  a  per  unit  of  surface  area  since 
the  other  heating  terms  in  the  dryer  model  are  represented  on  a  flux  basis  (i.e., 
heating  rate  per  unit  of  product  surface  area);  therefore, 


rxb 


PX 

AX 

o 


(4.20) 


A  representation  of  Equation  4.20  is  given  in  Figure  4.4A;  Qg  is  a  maximum  at  XQ 
and  then  decreases  linearly  to  zero.  Recent  experimentation  has  shown,  however, 
that  Qg  tends  to  be  constant  for  some  time  prior  to  this  linear  decline;  therefore. 
Equation  4.20  needs  to  be  modified  as  follows. 

If  the  product  moisture  level  is  higher  than  the  “saturation  point”  of  the 
dielectric  oven  (i.e.,  the  moisture  content  above  which  all  of  the  incident  power  is 
absorbed),  it  can  be  assumed  that  Qg  is  constant  The  hypothesis  in  Equation  4.18 
can  now  be  modified  to  reflect  this  observation.  Instead  of  using  the  ratio  Qv/X  , 
the  ratio  Qv/<5  is  used  where  OisX  normalized  by  the  “saturation  point”  (moisture 
content  at  which  the  drying  rate  first  begins  to  decline).  This  value  of  X  plays  the 
role  of  Xc  for  dielectric  drying  and  will  be  termed  Xcd.  Now,  the  volumetric  heat 
generation  per  unit  mass  of  water  is  constant  until  the  critical  moisture  content  is 
reached.  Following  the  same  procedure  as  before, 

Qv  =  r  =  p  =  (4.2i) 

o  OAb  oav 
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(4.22) 


For  <P  >  1,  Qg  is  constant;  for  <P  <  1  (i.e.,  below  Xcd),  Qg  =  T<I>.  Now  Qg  will 
remain  constant  until  Xcd,  and  then  will  fall  linearly  to  zero  for  X  <  Xcd,  as 
represented  in  Figure  4.4B. 

If  the  moisture  level  is  less  than  the  saturation  point  of  the  dielectric  oven, 
a  decay  of  the  absorbed  power  will  be  observed.  If  the  product  temperature  is  at 
equilibrium  (sensible  heat  is  negligible),  an  energy  balance  yields 
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where  the  primed  Qg  refers  to  the  dielectric  heat  flux  in  the  falling  rate  period. 
Therefore,  the  heat  generation  term  will  decay  proportionally  according  to 


Q'  =  —  Q 

g  (j)  g 


(4.25) 


In  the  dryer  model  it  is  proposed  that  Equation  4.25  be  used  for  the  falling  rate 
period  and  Equation  4.23  for  the  warm-up  and  constant  rate  periods  to  represent  the 
dielectric  heat  flux  as  a  function  of  the  material  properties. 


CHAPTER  5 
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THE  MODIFIED  DRYER  CODE 


The  dryer  code  developed  by  Flake  [2]  has  been  modified  to  handle 
dielectric  heating  and  boiling  point  drying.  The  basic  computational  procedure  used 
in  Flake's  model  is  essentially  unchanged.  It  is  convenient  to  take  X  as  the 
independent  variable  and  step  through  the  length  of  the  dryer  in  AX  increments 
from  the  known  initial  moisture  content  to  the  desired  exit  moisture  content, 
calculating  the  required  dryer  length  for  each  increment  The  original  dryer  code 
has  been  well  documented  in  Reference  [2];  thus,  only  significant  changes  made  in 
the  code  are  reported  here. 


The  most  important  subprograms  of  the  enhanced  code  are  described. 
The  code  is  validated  using  experimental  data  from  different  sources,  and  its 
performance  is  evaluated  by  means  of  running  hypothetical  cases  representative  of 
different  drying  regimes.  A  listing  of  the  code  and  a  variable  dictionary  are 
presented  in  Appendix  A. 
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Main  Program 


The  train  program  initializes  and  passes  all  the  variables  to  a  non-linear 
equation  solver,  called  NLSYST,  where  new  values  are  determined.  The  main 
program  writes  the  new  values  for  the  variables  and  re-initializes  them  for  another 
computational  step.  Figure  5.1  shows  the  relationship  between  the  main  program 
and  its  subroutines.  COMMON  blocks  were  extensively  used  to  avoid  transferring 
arguments  and  to  facilitate  the  “debugging”  of  the  Fortran  code. 


to 


This  subroutine  contains  all  the  governing  equations  describing  the  drying 
problem.  A  summary  of  these  equations  is  given  in  Appendix  B.  The  major 
modifications  made  to  this  subroutine  were: 

•  Incorporation  of  “IF’  statements  and  local  variables  that  reduced  the  sub-program 
to  half  of  its  original  length. 

•  Incorporation  of  a  subroutine  called  ASSIGN,  which  allows  governing  equations 
to  be  written  using  the  variables'  real  names  instead  of  meaningless  subscripted 
variables  as  required  by  the  subroutine  NLSYST. 


This  subroutine  represents  the  product  model  in  the  dryer  code.  When 
called,  it  computes  the  normalized  moisture  content,  the  relative  drying  rate,  and  the 
evaporative  front  position,  taking  into  consideration  whether  or  not  the  product  is  in 
the  warm-up  rate,  constant  rate,  or  falling  rate  period.  For  dielectrically  enhanced 
drying,  the  user  must  enter  a  flag  to  indicate  whether  the  behavior  of  the  drying 
curve  during  the  warm-up  rate  and  falling  rate  periods  is  linear,  parabolic,  or  cubic. 
These  functions  are  given  by  Equations  3.32-3.34,  and  are  already  incorporated  in 
the  code. 


This  program  models  the  dielectric  heating  within  the  wet  product. 
Equations  4.25  and  4.23  are  currently  used  in  the  code,  but  any  other  model  could 
be  used. 


MM 


5.3  VALIDATION  OF  THE  DRYER  CODE 
5.3.1  Introduction 

Validation  studies  of  the  original  dryer  code  were  carried  out  by  Rake  [2] 
for  purely  convective  cases.  To  validate  the  modified  code,  three  cases  were 
designed  in  which  experimental  data  from  the  literature  and  experimental  drying 
data  generated  in  ihe  Process  Energetics  Laboratory  of  the  University  of  Texas  at 
Austin  were  used  for  comparison  and  illustration.  The  cases  are  intended  to 
demonstrate  the  attributes  and  shortcomings  of  the  code.  Case  1  examines  the  code 
for  drying  without  an  evaporative  front.  The  experimental  data  for  this  case  was 
published  by  Lyons  [13]  et  al.  Case  2,  which  includes  an  evaporative  front,  gives 
a  step  by  step  explanation  in  how  to  apply  the  code  to  existing  experimental  data. 
Case  3  is  an  exercise  to  determine  the  amount  of  dielectric  power  necessary  with 
ambient  air  to  attain  the  same  drying  time  as  purely  convective  drying  with  air  at 
150°C. 


5.3.2  Casel:  Drying  without  an  evaporative  front 

Experimental  data  presented  by  Lyons  et  al.was  used  to  check  the 
modified  dryer  code.  The  material  used  was  absorbent  cotton  yam  with  the 
following  thermophysical  properties: 

Density,  ps  80  kg/m^ 

Thermal  conductivity,  Xs  0.00006  kW/m  K 

Specific  heat,  Cs  1.3kJ/kgK 

Surface  area,  A  0.0872  m^ 

Mass,  ms  1 .474  kg 

Air  conditions  were  not  documented;  Lyons  does  not  mention  an  air 
heater  when  discussing  his  experimental  apparatus.  He  attributes  minor  losses  to 
natural  convection  around  the  sample;  for  computational  purp r  -»s,  an  ambient  air 
velocity  of  0.4  m/s  was  specified.  The  ambient  air  temperaf  idity  and 

flow  rate  were  specified  to  avoid  condensation  of  moisture  ( 


Lyons  ct  al.  did  not  document  the  microwave  power  level;  therefore,  it 
was  estimated  from  the  temperature  history  of  the  sample.  This  approximation 
should  be  reasonably  accurate  at  the  beginning  of  the  drying  process  because  both 
the  air  and  the  wet  product  are  at  the  same  temperature,  thereby  minimizing 
convective  losses.  Deriving  3T/3t  from  the  temperature  history  graph  and  applying 
the  pertinent  conversion  factors,  the  temperature  rate  of  change  was 

=  0.0555  —  (5.2) 

3t  sec 

Performing  an  energy  balance  on  the  sample  and  neglecting  the  convective  term: 

3T 

P  =  Qw  A  =  Ms(C  +XCw)-^  =  0.8  kW  (5.3) 

This  value  is  consistent  with  the  size  of  the  commercial  microwave  oven  used  in  the 
tests. 

Data  normalization 

Under  ordinary  circumstances,  data  are  normalized  using  the  normalizing 
code,  NORM.  This  code  calculates  the  boiling  drying  rate,  Nv0,  the  drying  mode 
factor,  <)>,  and  the  relative  drying  rate,  /r .  The  relative  drying  rate  is  plotted  again?  ■ 
X  and  a  curve  fit  is  obtained.  Finally,  the  function  fT  vs.  X  is  input  into  the  dryer 
code  in  the  subroutine  "RMODEL".  A  different  approach  was  used  for  the  warm¬ 
up  period  in  this  case  study  because  the  temperature  history  of  the  product  as  a 
function  of  moisture  content  was  not  documented.  It  was  assumed  that  Nv  vs.  X 
was  linear  during  the  warm-up  period,  and  the  normalized  moisture  content  was 
defined  according  to  Equation  3.14. 


In  Lyons  et  al.'s  experiment  it  was  observed  that  sample  weight  loss  at 
the  beginning  of  the  experiment  was  due  to  water  “gushing”  out  and  not  to  actual 
evaporation.  This  “liquid  movement  period”,  as  they  defined  it,  is  not  modeled  by 


the  present  dryer  code;  thus,  it  was  decided  to  start  the  normalization  after  the  liquid 
movement  period.  The  X0  was  fixed,  and  Xb  was  iterated  on  until  Tp  =  100°C  and 
fx=  1  atX  =  XB. 


Results 


The  results  are  plotted  in  Figures  5.2  and  5.3.  The  predicted  drying  rate 
agreed  with  the  experimental  data  for  the  region  past  the  liquid  movement  period. 
The  experimental  data  for  the  liquid  movement  period  was  included  to  show  that 
even  though  the  simulation  was  started  at  a  lower  moisture  content,  the  slope  for 
the  drying  rate  curve  during  the  warm-up  period  was  similar  to  that  of  the 
experimental  drying  curve  as  it  was  established  in  Chapter  3. 

The  temperature  history  agreed  with  the  experimental  data  to  a  reasonable 
degree.  The  small  discrepancy  may  be  due  to  two  reasons:  first,  it  was  assumed  in 
the  code  that  the  whole  product  was  at  a  uniform  temperature;  and  second,  the 
experimental  temperature  data  plotted  in  Figure  5.2  was  taken  from  the  center  of  the 
product  which  may  be  different  from  an  average  temperature  across  the  sample. 

5.3.3  Case  2:  Drying  with  an  evaporative  front 

This  section  demonstrates  how  the  modified  code  is  applied  to  existing 
experimental  data  in  order  to  predict  drying  rates,  drying  times  and  dryer  length  for 
other  drying  conditions.  Experimental  data  were  collected  by  Armstrong  [24]  for  a 
polyurethane  foam  bed  with  a  surface  area  of  0. 1  m2  and  a  thickness  of  5  cm.  In 
this  test,  Ta  =150°C,  P  =  0.5k W,  and  Ua  =  3  m/s.  Temperature  and  sample  mass 
measurements  were  taken  as  explained  in  Section  3.3.  To  obtain  the  characteristic 
drying  function,  the  drying  rate  curve,  which  is  illustrated  in  Figure  5.4,  was 
normalized  according  to  the  procedure  listed  in  Section  3.3.3. 
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Figure  5.2  Temperature  history  prediction  for  Case  1.  Experimental  data  extracted 

from  Reference  [13]. 
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Figure  5.3  Drying  rate  prediction  for  Case  1.  Experimental  data  extracted  from 

Reference  [13]. 
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Several  key  points  were  selected  from  the  drying  rate  curve  and  the 
temperature  history  curve,  and  a  data  file  was  generated.  The  most  important 
points  to  consider  are  the  initial  and  final  conditions,  and  points  identifying  peaks 
and  curvatures  with  their  respective  surface  and  product  temperatures.  Presently, 
the  normalizing  code  looks  for  a  data  file  named  DATA,  in  which  the  values  of  X, 
Nv,  Ts  and  Tp  are  stored.  An  input  file,  which  is  presently  called  INPUT,  is 
generated  with  the  following  information  (see  variable  dictionary  in  Appendix  A  for 
definitions  and  units): 

A  B  POWER  N  EFC  TAIR 

VEL  NUC1  NUEXP1  LENGTH  KD  IPRINT 

For  the  case  study  the  following  conditions  were  used: 

0.1  0.05  0.5  300  0.05  150 

3  0.055  0.8  0.38  0.000021  10 

The  evaporative  front  coefficient  p,  (EFC)  was  set  to  0.05  since  capillary  forces 
are  expected  to  maintain  the  evaporative  plane  close  to  the  surface. 

Data  generated  by  the  normalizing  code  are  presented  in  Figure  5.5,  from 
which  Xq,  Xb,  Xcd,  f0,  and  /f,  the  data  points  necessary  to  develop  the  product 
model,  are  extracted.  The  characteristic  drying  curve  for  the  warm-up  and  falling 
rate  period  is  categorized  in  the  present  code  as  either  linear,  parabolic  or  cubic;  the 
shapes  of  these  functions  were  illustrated  in  Figure  3.6.  Based  on  Figure  5.5,  a 
linear  function  was  chosen  for  the  warm-up  period,  and  a  cubic  relation  was  chosen 
for  the  falling  rate  period;  the  cubic  function  (Equation  3.34)  is  plotted  with  the 
experimental  data  in  Figure  5.6. 

The  input  file  for  the  dryer  code  is  now  generated.  The  parameters  are 
read  in  the  following  order: 
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MWP 
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POWER 
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KD 
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For  the  Case  2  study,  the  input  parameters  were  determined  as  follows. 
Since  the  air  temperature  was  held  constant  at  150  °C,  ISO  was  set  to  1  to  represent 
isothermal  conditions  inside  the  dryer.  MWP  was  set  to  1  because  dielectric 
heating  was  used.  Delta  was  set  to  0.01.  SW  was  set  to  1  to  cause  the  product 
model  to  use  the  linear  characteristic  relation  between  X  and  Nv  for  the  warm-up 
period,  and  SF  was  set  to  3  to  cause  the  product  model  to  use  the  cubic 
characteristic  function  for  the  falling  rate  period  In  the  second  row,  G  was  set  to  a 
very  large  number  to  simulate  isohumid  conditions  inside  the  dryer.  L  was 
calculated  based  on  a  conveyor  belt  velocity  of  0.3  meters  per  minute.  The  third 
row  contains  the  parameters  associated  with  the  product  model.  These  parameters 
were  obtained  from  the  characteristic  drying  curve.  The  fourth  and  fifth  rows  are 
the  product  properties  and  convection  parameters,  respectively.  Documentation  of 
the  convective  parameters  is  given  by  Flake  [2],  Thus,  the  input  file  for  this  case 
is: 


1 

1 

300 

10 

0.01 

1 

3 

19.2 

0.06 

150 

25 

0.01 

30000 

0.01 

0.5 

10 

9 

15.54 

0.99 

0 

0.3 

0.2 

0.05 

3 

0.055 

0.8 

0.38 

0 

Note  that  the  units  used  in  the  code  are:  kg,  kJ,  kW,  °C,  seconds  and  meters. 

After  running  the  code,  an  output  file  is  generated  as  shown  in  Table  5.1. 
The  second  column  of  the  output  file  is  the  running  sum  of  the  required  dryer 
length  increments  as  the  moisture  content  is  decremented.  YA/YSAT  gives  the 
level  of  saturation  inside  the  dryer.  The  fifth  column  is  the  drying  mode  factor. 
Note  that  the  drying  process  is  dielectrically  dominated  towards  the  end.  This 
behavior  was  observed  in  all  the  runs  made. 

To  check  the  code  for  internal  consistency,  X  vs.  t  was  plotted  for  the 
experimental  data  and  the  dryer  model  data.  As  shown  in  Figure  5.7,  the 
agreement  is  very  close.  The  drying  time  estimated  by  the  code  was  within  about 
1  %  of  the  experimental  drying  time,  and  the  predicted  temperatures  were  also  in 


close  agreement  with  the  experimental  data,  as  illustrated  in  Figure  5.8. 
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DRYER  LENGTH  =  105.62  METERS 
GENERATED  POWER  =  5.00  KW/SQ.METER 


DRYING  TIME  =  337.97  MINUTES 
AIR  TEMP.=  150.0  DEGREES  C 


Table  5. 1  Output  file  for  Case  2 
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Figure  5.7  Moisture  content  history  for  Case  2. 
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Figure  5.8  Temperature  history  as  a  function  of  moisture  content  for  Case  2. 


Case  3  considers  the  same  drying  conditions  and  material  in  Case  2,  but 
without  dielectric  heating.  A  drying  test  was  performed  by  Armstrong  [24]  for 
these  conditions;  it  took  approximately  57  hours  for  the  foam  to  dry  from  2000%  to 
0  %  moisture..  The  objective  of  this  case  is  to  estimate  the  dielectric  power  that 
would  be  necessary  to  achieve  the  same  drying  time  using  ambient  temperature  air. 
Ta  is  set  to  25°C,  MWP  is  set  to  one,  and  the  dielectric  power  is  iterated  until  the 
same  drying  time  is  attained.  For  this  case  the  code  predicted  that  the  required 
power  would  be  1.47  kW/m2. 

This  case  illustrates  one  of  the  advantages  of  dielectric  heating  for  drying. 
For  the  convective  case,  the  wet  product  reached  a  maximum  temperature  of  70°C, 
while  the  surface  temperature  reached  125  °C.  For  the  hybrid  case,  the  dryer  code 
predicted  that  the  wet  product  will  reach  a  maximum  of  100°C  while  the  surface 
temperature  only  reaches  73.5°C.  Dielectric  heating  can  play  an  important  role  in 
drying  temperature  sensitive  materials,  producing  equivalent  drying  rates  at  lower 
temperatures  than  pure  convective  drying. 

Drying  can  also  be  accelerated  using  higher  power  densities,  still  without 
requiring  high  air  temperature.  To  see  the  effect  of  dielectric  heating  on  the  drying 
time,  the  dielectric  power  input  was  doubled  while  maintaining  Ta  =  25°C.  The 
drying  time  was  reduced  to  818  minutes,  a  76%  reduction,  but  the  predicted  surface 
temperature  increased  by  only  0.5°C. 


5.4  DRYER  PERFORMAN 


HEBE 


ER  VARYING  DRYING  REGIME 


Several  hypothetical  cases  were  specified  based  on  the  drying  regimes 
defined  in  section  3.2  to  test  the  code  under  different  drying  conditions.  Table  5.2 
is  a  compilation  of  the  results.  The  first  three  columns  in  the  tabie  are  the  initial 
drying  conditions  relevant  to  dielectric  heating.  The  maximum  surface  temperature, 


maximum  wet  product  temperature,  and  drying  time  are  tabulated  for  comparison 
purposes.  The  first  row  of  Table  5.2  (pure  convective)  is  experimental  data 
published  by  Armstrong  [24]  as  discussed  in  Case  3  above. 
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The  trends  observed  in  the  drying  times  and  temperatures  behaved  as 
expected.  The  most  interesting  result  occurs  in  the  subcooled  dielectrically 
dominated  region.  The  predicted  drying  time  was  more  than  twenty  times  lower 
than  the  drying  time  for  the  purely  convective  case  and  the  convectively  dominated 
case;  however,  the  surface  temperature  reached  a  temperature  of  78°C,  considerably 
lower  than  the  other  two  cases. 

Results  could  not  be  obtained  for  the  subcooled  convectively  dominated 
case  due  to  instability  problems  with  the  numerical  nonlinear  equation  solver.  For 
the  pure  convective  cooling  regime  and  the  dielectrically  dominated  cooling  regime 
the  code  behaved  erratically.  The  number  of  computational  steps  was  increased 
with  no  positive  results.  To  simulate  the  cooling  regime,  the  initial  product 
temperature  was  set  to  80°C,  and  while  the  code  did  run  for  these  conditions,  it 
generated  unreasonable  data,  indicating  instability  for  these  drying  conditions. 

5.4.2  Vapor  flux  effect  in  drv  region 

The  code  was  modified  to  account  for  the  effect  of  vapor  flux  on  the  dry 
region  temperature  profile  as  discussed  in  Section  4.1.2.  To  determine  the  effect  of 
this  change,  the  code  was  run  with  and  without  the  modification  using  the  drying 
conditions  specified  for  Case  2,  but  at  a  much  higher  microwave  power  level.  The 
results  are  tabulated  in  the  last  two  rows  of  Table  5.2.  There  was  a  small  difference 
in  the  maximum  predicted  surface  temperature  and  a  large  difference  in  the  final 
surface  temperature  at  the  end  of  the  drying  process.  Despite  these  differences,  the 
drying  time  was  not  significantly  affected.  This  is  not  surprising,  since  at  high 
dielectric  heating  levels,  the  convective  cooling  term,  which  is  the  one  most  affected 
by  the  predicted  product  surface  temperature,  is  negligible. 
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REGIME 
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CHAPTER  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  and  recommendations  presented  in  this  section  are 

associated  with  development  of  predictive  tools  for  evaluation  of  applications  of 

dielectrically  enhanced  dryers. 

6.1  CONCLUSIONS 

In  general,  it  is  concluded  that  the  characteristic  drying  curve  concept 

should  be  applicable  to  dielectric  drying.  Specifically  in  this  study: 

•  A  unique  characteristic  drying  curve  was  obtained  for  dielectrically  assisted 
drying  of  a  bed  of  glass  beads. 

•  A  dimensionless  energy  balance  equation  was  formulated  for  dielectrically 
assisted  drying,  and  used  to  define  a  dimensionless  parameter,  the  drying  mode 
factor,  which  was  used  to  define  drying  regimes. 

•  Relative  drying  rate  and  relative  warm-up  rate  functions  were  defined  which 
depend  only  on  the  moisture  content  and  the  material  properties. 

•  An  equation  was  developed  to  estimate  the  position  of  the  evaporative  front  in 
drying  capillary-porous  materials. 

•  Governing  equations  were  modified  to  account  for  the  thermal  effect  of  vapor 
transport  in  the  product's  dry  region. 

•  A  heat  generation  term  model  was  developed  and  incorporated  into  the  dryer 
code. 
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•  The  existing  computer  code  was  modified  to  accept  dielectrically  assisted  drying. 
Sanity  checks  were  performed  on  the  computer  code  by  comparing  the  results 
obtained  using  the  code  to  published  data  found  in  the  literature  and  to 
experimental  data  attained  in  the  Process  Energetics  Laboratory  at  the  University 
of  Texas  at  Austin. 


It  was  concluded  from  the  dimensional  analysis  of  the  energy  balance 
equation  that  the  relative  drying  rate  and  the  relative  warm-up  rate  were  independent 
of  power  level  and  air  temperature.  This  concept  was  experimentally  confirmed  by 
applying  the  normalizing  scheme  to  drying  data  taken  at  two  microwave  power 
levels  (0.5  kW,  1.0  kW)  with  different  air  temperatures  (65°C,  100 °C,  120°C, 
150°C). 

6.2  RECOMMENDATIONS 

Presently,  the  code  uses  the  conditions  at  the  entrance  of  the  dryer  as 
input;  for  countercurrent  air  flow  the  air  conditions  are  unknown  at  the  dryer 
entrance,  and  therefore,  the  exiting  air  temperature  must  be  estimated  until  the 
computed  initial  air  conditions  match  the  given  initial  air  conditions.  Currently  this 
iteration  is  performed  manually.  It  is  recommended  that  the  necessary 
programming  instructions  be  incorporated  into  the  code  so  that  iterations  are 
performed  automatically. 

It  is  important  to  further  modify  the  dielectrically  enhanced  drying  model  to 
handle  multi-zone  drying.  Drying  data  are  needed  to  evaluate  the  behavior  of  the 
drying  rate  in  the  transition  from  dielectric  to  conventional  drying  and  vice  versa. 

In  principle,  each  drying  zone  can  be  treated  as  a  different  drying  problem  whose 
input  data  is  obtained  from  the  exit  conditions  of  the  prior  zone.  The  main 
difficulty  rests  in  developing  a  product  model  to  handle  the  transition  from  one  zone 
to  the  other;  in  other  words,  the  product  model  must  be  independent  of  the  moisture 
content  to  allow  the  code  to  compute  a  relative  drying  rate  from  the  characteristic 
drying  function  assigned  to  a  particular  drying  zone.  The  enhanced  code  developed 


here  has  been  structured  in  modular  form  to  permit  extension  to  the  multi- zone 
problem.  This  assumes,  however,  that  the  material’s  characteristic  drying  curve 
will  respond  instantaneously  to  a  sudden  change  in  boundary  conditions.  This 
question  has  not  been  investigated. 

The  characteristic  drying  curve  also  needs  to  be  correlated  to  the  initial 
moisture  content  and  initial  product  temperature.  These  two  parameters  are  very 
significant  for  the  development  of  a  multi-zone  dryer  model  because  the  initial 
conditions  are  going  to  be  different  for  each  different  zone  encountered  by  the 
product. 

It  is  recommended  that  <)>  be  used  as  a  management  tool  for  the  planning  of 
dielectric  drying  experiments.  Existing  hardware  would  have  to  be  modified  to 
allow  lower  microwave  power  levels  and  higher  air  temperatures  in  order  to  achieve 
the  convection  dominated  regime  (0  >  2)  and  the  subcooled  regime  (0  <  1). 
Knowing  more  about  the  convective  regime  is  important  because  it  concerns  the 
incorporation  of  dielectric  power  into  existing  convective  dryers.  The  subcooled 
regime  is  of  interest  because  it  permits  drying  with  ambient  air.  Table  6.1  shows  a 
proposed  matrix  for  the  experimental  runs.  AT  is  the  equilibrium  temperature 
difference  between  the  air  and  the  product  surface.  The  table  is  symmetrical  about 
the  parameter  a,  the  ratio  of  Qfc/Qg,  where 

0  =  1  +  a  (6.1) 

Note  that  a  is  negative  for  convective  cooling  and  positive  for  convective  heating 
which  generates  a  mirror  image  of  the  table  to  the  left  of  AT  =  0.  Some  different 
combinations  of  AT  and  P  provide  the  same  0,  so  while  Table  6.1  reflects  55 
experiments,  only  37  are  necessary  to  cover  the  range  of  0's.  Values  for  P  and  AT 
should  be  selected  according  to  the  experimental  equipment  capacity.  To  generate 
Table  6.1,  it  was  assumed  that  the  convective  heat  transfer  coefficient  stays 
relatively  constant  when  going  from  -5AT  to  5AT. 
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1±  3a 

l±4a 

l±5a 

2P 

1 

1±  0.5a 

1±  a 

1+  1.5a 

1±  2a 

1±  2.5a 

3P 

1 

1±  0.33a 

1±  0.67a 

1±  a 

1±  1.33a 

1±  1.67a 

4P 

1 

1±  0.25a 

1±  0.5a 

1±  0.75a 

1±  a 

1±  1.25a 

5P 

1 

1±  0.2a 

1±  0.4a 

1±  0.6a 

1±  0.8a 

1+  a 

Table  6. 1  Drying  mode  factor  matrix  for  dielectrically  enhanced  drying  tests. 


The  present  code  provides  critical  input  information  for  economic 
comparison  of  conventional  and  hybrid  dryers.  All  parameters  that  might  affect  the 
ultimate  economic  efficiency  ($/Kg  of  removed  water)  of  the  dryer  should  be  taken 
into  account.  As  mentioned  in  Chapter  2,  Pourhiet  [12]  developed  an  expression  to 
convert  dielectric  drying  transfer  units  in  a  hybrid  system  to  conventional  drying 
transfer  units.  A  similar  concept  could  be  incoiporated  into  this  model  to  aid  in 
comparing  the  two  approaches. 

One  of  the  assumptions  made  when  developing  the  model  was  that  a 
linear  temperature  profile  exists  in  the  dry  region.  This  implies  that  the  dry  region’s 
thermal  capacitance  is  small.  This  assumption  is  not  valid  when  the  thermal 
capacitance  of  the  product  is  relatively  high.  Not  knowing  the  temperature  profile 
inside  the  dry  region  precludes  estimating  the  amount  of  heat  being  conducted 
through  the  dry  region  into  the  wet  region.  It  is  recommended  that  the  dry  region 
model  be  modified  by  discretizing  it  in  the  y  direction  and  by  using  the  heat 
diffusion  equation  to  determine  the  temperature  profile,  subject  to  the  following 
boundary  conditions: 


T  =  T]  , 
p  Jy=5 


5y  Jy=0 


After  discretizing  the  diffusion  equation,  n  equations  in  n  unknowns  are  added  to 
the  already  existing  equations  in  the  dryer  model.  The  calculation  of  the 
temperature  profile  inside  the  dry  region  could  be  incorporated  into  an  optional 
subroutine.  An  alternative  approach  has  been  suggested  in  which  a  comprehensive 
physical  predictive  product  model  is  interfaced  with  the  empirical  approach  to 
generate  a  characteristic  drying  curve  for  convective  and  enhanced  drying.  Such 
models  are  being  developed  as  part  of  the  overall  drying  research  program  at  the 
University  of  Texas  at  Austin. 


The  model  used  to  predict  the  position  of  the  evaporative  front  must  be 
experimentally  validated.  Knowing  8  is  crucial,  especially  when  drying  thick  low 
thermal  conductivity  materials  with  air  at  temperatures  much  higher  or  lower  than 
the  boiling  temperature. 

The  validation  studies  performed  on  the  dryer  code  produced  positive 
results  within  the  limited  range  of  available  drying  data.  Experiments  need  to  be 
run  to  test  the  predictive  capability  of  the  modified  code  for  a  wider  range  of  drying 
conditions.  For  example,  a  drying  test  should  be  performed  on  polyurethane  foam 
for  the  same  initial  moisture  content  and  initial  product  temperature  as  in  Case  2, 
but  at  a  microwave  power  level  of  15  kW/m2  and  air  temperature  of  65°C.  A 
characteristic  drying  function  for  foam  is  already  documented  in  this  report;  an 
immediate  follow-on  should  be  to  confirm  its  ability  to  extrapolate  to  conditions 
beyond  the  range  of  the  data  used  to  generate  the  function. 

The  problem  of  numerical  instability  of  the  code  during  cooling  of  the 
product  must  be  investigated  and  resolved.  It  is  unclear  whether  the  problem  lies  in 
the  inadequacy  of  the  nonlinear  equation  solver  routine  or  in  inconsistency  of  the 
basic  product  model  in  the  cooling  regime.  Since  cooling  conditions  might 
commonly  occur  in  multi-zone  dryers,  this  question  is  critical  to  successful 
development  of  a  multi- zone  model. 
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A.l  VARIABLES  DICTIONARY 


A.2  LISTING  OF  THE  ENHANCED  DRYER  CODE 


A.3  LISTING  OF  THE  NORMALIZING  CODE 


r-- 

to 


K< 


A 

ASSIGN 

ATM 

B 

BETA 

CPA 

CPAF 

CPAI 

CPF 

CPFA 

CPFF 

CPS 

CPV 

CPVA 

CPVF 

D 

DAB 


-  Products  surface  area,  (m2) 

-  Subroutine  that  assigns  a  subscripted  variable  to  the  unknowns  as 
required  by  NLSYST  and  viceversa. 

-  Air  pressure  in  dryer,  (atmospheres). 

-  Product  thickness,  (meters) 

-  Coefficient  used  in  heat  and  mass  transfer  analogy,  (dimensionless) 

-  Dry  air  specific  heat,  (kJ/kg°C) 

-  Function  subroutine  used  to  calculate  dry  air  specific  heat,  CPA,  as 
a  function  of  temperature. 

-  Dry  air  specific  heat  under  isothermal  drying  conditions.  This  is 
separately  defined  because  it  remains  constant  throughout  the  dryer. 

-  Specific  heat  of  saturated  liquid  water,  (kJ/kg°C) 

-  Same  as  above,  evaluated  at  Tdp  for  use  in  air  enthalpy  equation, 

-  Function  subroutine  used  to  calculate  saturated  liquid  water  specific 
neat 

-  Specific  heat  of  the  (solid)  product,  an  input  variable,  (kJ/kg°C). 

-  Specific  heat  of  water  vapor,  (kJ/kg  K) 

-  Same  as  above,  evaluated  at  Ta  for  use  in  air  enthalpy  equation, 

-  Function  subroutine  used  to  calculate  water  vapor  specific  heat 

-  Ratio  of  molecular  weights  of  air  to  water,  (dimensionless). 

-  Diffusivity  of  water  vapor  in  air,  Dab,  (m2/s). 

-  Same  as  above,  except  computed  according  to  Eq.  B-12  of 
Reference  [2], 


DAB1 


DELTA  -  Input  variable  used  to  calculate  perturbation  value  in  nonlinear 

equation  solver  subroutine,  NLSYST,  where  for  a  given  unknown 
variable  X(I),  the  perturbation  amount  for  that  variable  is 
DELTA*X(I).  See  discussion  at  end  of  Chapter  3  of  Reference  [2] 

DELTAX  -  Equivalent  to  DELTX  (below),  except  at  end  nodes  DELTAX  = 
DELTX/2. 

DELTX  -  Decrement  value  of  moisture  content,  (AX),  equivalent  to  Xjn- 
Xout/N,  where  N  is  the  number  of  nodes,  (kg/kg). 

DELTY  -  Air  humidity  increment  (or  decrement  in  countercurrent  case), 

...  ere  DELTY  =  (L/G)  DELTX,  (kg/kg). 

DELZ  -  Increment  value  of  dryer  length,  Az,(m)  See  ZPOS. 

DMF  -  Drying  mode  factor,  <j>,  (dimensionless) 

DRYTTM  -  Drying  time,  function  of  dryer  position  and  product  flow  rate,  (s). 

EFC  -  Evaporative  front  coefficient,  P  (dimensionless) 

EFG  -  Change  of  5  with  respect  to  X,  (m-kg  solid/kg  water) 

ENTHAL  -  Humid  air  enthalpy,  ha,at  i*  step,  (kJ/kg). 

ENTHO  -  Humid  air  enthalpy  at  i-1  step,  set  equal  to  computed  air  enthalpy, 

ENTHAL,  after  each  step,  (kJ/kg). 

EQUELF  -  Equilibrium  moisture  content  factor.  Used  to  determine  the 

equilibrium  moisture  content  for  hygroscopic  materials,  which  is  the 
product  of  the  equilibrium  factor  and  the  air  humidity. 

EXPON  -  Exponent  n,  in  the  function  for  the  characteristic  drying  curve,  /  = 

On. 

F  -  Array  of  equations  which,  with  X(I)  unknowns,  is  solved  by 

NLSYST. 

F(l)  -  Mass  balance,  composite  CV  I,  CV II. 

F(2)  -  Energy  balance,  air  control  volume,  CV  I. 

F(3)  -  Air  enthalpy  equation. 


FCN 

FF 

FI 

FKAIR 

FRACT 

FREL 

FT 

FTOL 

G 

HC 


F(4)  -  Dew  point  temperature  function. 

F(5)  -  Energy  balance  for  adiabatic  saturation  process. 

F(6)  -  Saturated  humidity  function. 

F(7)  -  Energy  balance,  wet  region,  CV  HI. 

F(8)  -  Energy  balance,  dry  region,  CV  II. . 

F(9)  -  Drying  rate  equation. 

F(10)  -  Surface  energy  balance  for  wet  bulb  conditions. 

-  Subroutine  containing  equations  of  heat  and  mass  balances,  etc, 
solved  at  each  step  by  NLSYST. 

-  Final  relative  drying  rate.(dimensionless) 

-  Initial  relative  drying  rate.(dimensionless). 

-  Function  used  to  compute  thermal  conductivity  of  air,  ka,  (kW/m 
K). 

-  Relative  depth  of  evaporative  plane  (dimensionless). 

-  Relative  drying  rate  /r. (dimensionless) 

-  Array  of  dew  point  and  wet  bulb  equations  solved  by  initial  call  to 
NLSYST. 

-  Tolerance  parameter  used  by  NLSYST  to  determine  convergence. 
As  a  solution  converges,  the  equations  being  solved  F(I)  approach 
zero.  Convergence  is  reached  when  F(I)  <  FTOL.  See  Reference 
[2]. 

-  Air  mass  flow  rate,  dry  air  weight  basis,  (kgdry  au/s-m).  Positive 
value  indicates  concurrent  air  flow,  and  a  negative  value, 
countercurrent. 


HCF 


-  Convective  heat  transfer  coefficient,  h  (kW/m2°C) 

-  Function  subroutine  used  to  calculated  convective  heat  transfer 
coefficient 


HFGA 


-  Latent  heat  of  vaporization  evaluated  at  T<jp,  used  in  air  enthalpy 
equations. 


HFGS  -  Same  as  above,  evaluated  at  wet  bulb  temperature  used  in  wet  bulb 
equation. 

HFGW  -  Same  as  above,  evaluated  at  adiabatic  saturation  temperature,  Tas, 
used  in  adiabatic  saturation  energy  balance  equation. 

HFGF  -  Function  used  to  calculate  latent  heat  of  vaporization  of  water. 

I  -  Counter,  steps  from  1  to  N+l. 

IFLAG  -  Flag  set  when  air  humidity  is  greater  than  saturation  value. 

Saturation  will  occur  when  air  flow  rate  is  not  sufficient  to  entrain 
moisture  removed  from  product,  when  specified  entering  air 
humidity  is  at  or  greater  than  saturated  humidity  or  when  entering 
product  temperature  is  less  than  the  dew  point  temperature 
(condensation). 

IP  -  Printing  counter  reset  to  1  after  each  print.  See  IPRINT  below. 

IPRINT  -  Counter  used  to  determine  printing  increments.  For  example,  if  IP 

is  set  to  4,  variable  are  printed  every  fourth  step  (1,5,9,13, 

. ,N+1). 

ISO  -  Flag  set  to  1  if  program  is  to  be  run  in  an  isothermal  configuration, 

set  to  2  if  drying  configuration  is  nonisothermal 

K  -  Flag  sent  to  NLSYST.  If  set  to  one,  values  of  the  solved  variables 

and  equations  will  be  printed  at  each  iteration.  If  solution  fails  to 
converge,  the  equations  or  variables  which  did  not  converge  may  be 
found  by  examining  the  output.  As  currently  configured,  these 
values  are  printed  only  on  the  first  call  to  NLSYST. 

KD  -  Thermal  conductivity  of  dry  product,  an  input  variable,  (kW/m°C) 

L  -  Mass  flow  rate  of  product,  dry  mass  basis,  (kgdry  product/s-m) 

LE  -  Lewis  number. 

LENGTH  -  Characteristic  length  used  in  Nusselt  number  correlation,  an  input 
variable. 

MAXIT  -  Maximum  number  (limit)  of  iterations  to  be  made  in  NLSYST. 

MWP  -  Flag  set  to  1  if  dielectric  power  is  on;  set  to  zero  for  conventional 

drying. 


N  -  Number  of  nodes,  input  variable. 

NEQ  -  Number  of  equations  to  be  solved  by  NLSYST 

NTQ  -  Number  of  equations  to  be  solved  by  NLSYST  in  initial  call  for 

finding  wet  bulb  and  dew  point  temperatures  (NTQ=2). 

NTU  -  Number  of  transfer  units. 

NUC1  -  Coefficient ,  C.used  in  Nusselt  number  correlation. 

NUEXP1  -  Exponent ,  m,  of  Reynolds  number  used  in  Nusselt  number 

NLSYST  -  Subroutine  used  to  solve  a  system  of  nonlinear 
equations. 

NUL  -  Nusselt  number. 

NV  -  Drying  rate  flux,  Nv,(kg/m2s). 

NVO  -  Maximun  drying  rate  if  all  the  dielectric  heat  is  used  for  evaporation. 

Convective  heat  is  not  included,  (kg/sec-m2) 

PA  -  Air  pressure  in  Pascals. 

PD  -  Dielectric  power  deposition  in  the  wet  material,  (kW/m3). 

PEN  -  Argument  used  in  subroutine  WRITER.  When  set  to  zero,  the 

headings  of  the  output  files  are  writen;  when  set  to  one,  the 
variables  for  each  moisture  content  decrement  are  writen;  when  set 
to  two,  the  dryer  length,  drying  time  and  NTU  are  written. 

PHIH  -  Correction  to  convective  heat  transfer  coefficients  due  to  high  mass 
transfer  rates. 

POWER  -  Forward  dielectric  power,  (kW) 

PR  -  Prandtl  number. 

PVEL  -  Velocity  of  product  through  dryer,  (m/s). 

PWA  -  Partial  pressure  of  water  vapor  in  air,  (Pascals). 

-  Saturated  partial  pressure  of  water  vapor  (Pascals). 


PWS 


79 


QAS  -  Air  heating  source  flux  Qas»(kW/m2),  an  input  variable.  For 

isothermal  configuration,  the  air  heating  necessary  to  provide 
constant  temperature  air  is  calculated  and  the  input  value  of  QAS  is 
not  used. 

QCOND  -  The  heat  conducted  through  the  dry  region  (kW/m2). 

QCONV  -  The  heat  convected  at  the  surface  (kW/m2). 

QEVAP  -  The  total  rate  of  latent  heat  of  vaporization,  Nv  Ahv  (kW/m2). . 

QGEN  -  Dielectric  heat  flux,  (kW/m2) 

QGENF  -  Subroutine  used  to  calculate  QGEN 

QSENSD  -  Sensible  heating  in  the  dry  region  control  volume,  (kW/m2). 

QSENSW  -  Sensible  heating  in  the  wet  region  control  volume,  (kW/m2). 

QSOL  -  Enthalpy  flux  carried  in  by  the  dry  solid  due  to  the  displacement  of 
the  evaporative  front,  (kW/m2) 

QVAP  -  Change  in  enthalpy  of  the  water  vapor  when  transported  through  the 
dry  region,  (kW/m2) 

REL  -  Reynolds  number,  ReL.  used  in  Nusselt  number  correlation. 

RELHUM  -  Relative  humidity,  (dimensionless). 

RHOA  -  Air  density,  ps,  (kg/m3). 

RHODS  -  Density  of  the  dry  product  (solid),  an  input  variable,  (kg/m3). 

RKO  -  The  mass  transfer  coefficient,  Kq  ,  (Kg/m2s) 

RKOF  -  Function  subroutine  used  to  calculate  mass  transfer  coefficient  IQ,. 

RMODEL  -  Product  model  subroutine  used  to  calculate  relative  drying  rate,  /r, 
as  a  function  of  normalized  moisture  content. 

SC 


SW 


-  Schmidt  number. 

-  Flag  used  to  specify  whether  the  characteristic  function  during  the 
warm-up  period  is  linear,  parabolic  or  cubic. 


SF  -  Rag  used  to  specify  whether  the  characteristic  function  during  the 

falling  rate  period  is  linear,  parabolic  or  cubic. 

TA  -  Air  temperature,  Ta,  in  nonisothermal  configuration  (°C). 

TADSAT  -  Adiabatic  saturation  temperature,  Tas,  (°C). 

TAI  -  Air  temperature  in  isothermal  case,  constant  throughout  the  dryer, 

(°C). 

TAIR  -  Local  variable  used  in  FCN  to  represent  air  temperature,  (°C). 

TAL  -  Air  temperature  at  product  inlet  Qeft),  an  input  variable,  (°C).  This 

is  the  entering  air  temperature  in  concurrent  configuration,  exiting 
air  temperature  in  countercurrent. 

TAT  F  -  Same  as  above  but  in  °F.  Used  in  initial  calculations  of  the  wet  bulb 
temperature. 

TDEW  -  Dew  point  temperature,  Tdp,  (°C). 

TEVAPO  -  Evaporative  plane  temperature  at  previous  iteration  (old),  (°C) 

TFCN  -  Subroutine  containing  functions  for  dew  point  and  wet  bulb 
temperature,  called  by  NLSYST. 

TGU  -  Initial  guess  of  wet  bulb  temperature  (°F),  before  calling  NLSYST 

to  solve  for  wet  bulb  and  dew  point  temperatures. 

TGUC  -  Same  as  above,  but  in  °C. 

TMEANO  -  Mean  temperature  in  dry  product  region  (°C),  at  previous  iteration 
(old). 

TNITIAL  -  Subroutine  that  sets  flags  and  initializes  variables. 

TP  -  Product  temperature,  Tp,  (°C). 

TPL  -  Product  temperature  at  product  inlet  (left),  an  input  variable  (°C). 

TPOLD  -  Product  temperature  at  previous  iteration,  (°C). 

TSURF  -  Product  surface  temperature,  Ts,  (°C). 

TSURFO  -  Product  surface  temperature  at  previous  iteration,  (°C). 


TWB 


-  Wet  bulb  temperature,  Twb,  (°C) 


TWMIN 

VEL 

vise 

VISCF 

WETBULB 

WRITER 

X 

XB 

XCD 

XCR 

XMC 

XMCL 

XMCR 

XNORM 

XT 

XTOL 

YA 

YAL 

YS 


Minimum  wet  bulb  temperature  (°F),  used  in  initial  guess  calculation 
for  wet  bulb  temperature. 

-  Air  velocity,  an  input  variable  (m/s).  Used  to  compute  Reynolds 
number  in  heat  transfer  correlation  calculations. 

-  Air  viscosity,  kinematic,  va,  (m2/s). 

-  Function  used  to  calculated  air  viscosity,  defined  above. 

-  Subroutine  that  estimates  the  wet  bulb  and  adiabatic  saturation 
temperature. 

-  Subroutine  containing  write  and  format  statements  for  output  files. 

-  Variable  array,  X(I),  representing  the  unknown  variables  to  be 
solved  with  equations  F(I),  given  in  subroutine  FCN. 

-  Moisture  content  at  which  boiling  point  is  reached. 

-  Critical  moisture  content,  XC(j,  for  dielectric  heating;  an  input 
variable,  (kg/kg). 

-  Critical  moisture  content,  Xc,  an  input  variable,  (kg/kg). 

-  Product  moisture  content  at  current  step,  X,  (kg/kg). 

-  Product  moisture  content  at  inlet  (left),  an  input  variable  (kg/kg). 

-  Desired  product  moisture  content  at  outlet  (right),  an  input  variable 
(kg/kg). 

-  Normalized  moisture  content,  C>. 

-  Variable  array  where  XT(1)  is  Tdp,  and  XT(2)  is  Twb-  Solved  with 
equation  array  FT  by  NLSYST. 

-  Tolerance  value  used  by  NLSYST  to  determine  convergence. 
Iteration  is  stopped  at  X(I)-X(I)0i(j<XT 

-  Current  air  humidity,  Y,  (kg/kg). 

-  Air  humidity  at  product  inlet  (left),  an  input  variable,  (kg/kg). 

-  Surface  humidity,  Ys,  (kg/kg). 

-  Saturated  air  humidity,  Ysat,  (kg/kg). 


YSAT 


YWF 

YWSAT 


-  Function  subroutine  used  to  calculated  saturated  humidity  as  a 
function  of  temperature. 

-  Saturated  air  humidity  evaluated  at  wet  bulb  temperature,  YW) 
(kg/kg). 


ZPOS 


-  Current  position  in  dryer,  measured  from  the  product  inlet,  (m). 
Computed  by  summation  of  A z. 
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A. 2  LISTING  OF  THE  ENHANCED  DRYER  CODE 


PROGRAM  DRYER  (INPUT, TAPE4,TAPE6,TAPE7,TAPE8,TAPE5,TAPE9=  INPUT, 

*  TAPE10,TAPE11,TAPE12,TAPE1,TAPE2,TAPE3,TAPE13,TAPE14,TAPE15, 

*  TAPE  1 6.TAPE  1 7,T APE  1 8, TAPE  1 9.TAPE20.T APE2 1  .TAP E22.T APE23 . 

*  OUTPUT, TAPE24=OUTPUD 

I!************************************************************************* 


ENHANCED  TUNNEL  DRYER 

THIS  PROGRAM  SIMULATES  PROCESS  CONDITIONS  IN  A 

CONTINUOUS  CONVECTIVE,  DIELECTRICALLY  ENHANCED  TUNNEL  DRYER 

UNDER  ADIABATIC  OR  ISOTHERMAL  AIR  CONDITIONS. 

AT  PRESENT,  THE  CODE  SUCCESSFULLY  HANDLES  PRODUCT 
MODELS  OF  PURELY  CONVECTIVE  OR  HIBRIDLY  DRIED  MATERIALS  AND  IS 
STRUCTURED  SUCH  THAT  FUTURE  EDITIONS  CAN  INCORPORATE 
EQUIPMENT  MODELS  FOR  MULTI-ZONE  DRYING. 

THE  ACCOMPANYING  VARIABLE  DICTIONARY  PROVIDES  INFORMATION 
ON  VARIABLES  AND  SUBROUTINES  FOR  USE  OF  THE  CODE. 

************************************************************************* 


EXTERNAL  FCN.TFCN 
DIMENSION  F(16),X(16),XT(2),FT(2) 

REAL  KO.LLELENGTHNUCINUEXPINULNTUNVNVO 

COMMON/FUNCT/DELZ.DR  YTIM  £NTHAL£NTHOOMF.NEQ,NTU  JMV.NVO.PD, 

*  QAS,QCOND,QCONV,QEVAP,POWER,QGEN,QSENSD,QSENSW,TA,TADSAT, 

*  TDEW.TMEANO.TP,TPOLD,TSURF.TSURFO,TWB,YSAT.ZPOS.QVAP,QSOL 
COMMON  A.B,D,DELTAX,GJJSO,LMWP,PWS,RELHUM,TAIXMC,PA,YA 
COMMON/NLSYS/MAXITX)ELTAXrOLFTOL 

COMMON/PROP/CP  A,CPAI,CPF,CPFACPS,CPV,CPVA,FKA,HFGA,HFGS, 

*  HFGW.KDJIHODS.VISC 

COMMON/PMODElVEQUILF3W,EXPONJFRACT,FRELXCRANORM,SF,EFC,EFG, 

*  FLXB.XCD,FF 

COMMON/CONV/ATM3ETA,CPRAT,CP  Y.DAB  ,HCLE,LENGTUNUC  1  ,NUEXP1,NUL 

*  PHIH,PR,PSYRAT,RELRHOA,RHOMA,RKO,SC,VELYDA,YDS,YS 
COMMON/TTNIT/IFLAG 1  JFLAGZK,YALXMCLXMCR,TPLTALJFLAG 

CALL  WRITER(0.) 

D  =  0.62198 

- READ  INPUT  PARAMETER. STATEMENTS  ARE  DIVIDED  AS — 

(1) FLAGS, 

(2) INmAL  CONDITIONS,  DRYER  CONDITIONS 

(3)  PRODUCT  MODEL'S  COEFFICIENTS, 

(4)  THERMOPHISICAL  PROPERTIES, 

(5)  CONVECTIVE  PARAMETERS. 

READ  (9,*)  ISO.M  AT.NJPRINT.DELTA.SW'F 


non  noon  non 
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READ  (9,*)  XMCL,XMCR,TAL,TPL,YAJUG,L,  POWER 
READ  (9,*)  XCR.XCD,XB,EXPON,EQUILF,FI,FF,EFC 
READ  (9,*)  CPS.PA.B.KD.RHODS  A 
READ  (9,*)  VEUNUC1.NUEXP1, LENGTH, QAS 
C 

C - CALCULATE  AX  AND  AY,  PRESSURE  IN  PASCALS - 

C 

IP=  IPRINT 

DELTX  =  (XMCR-XMCLyFLOAT(N) 

DELTY  =  -(LG)*  DELTX 
DELTAX  =  DELTX 
PA  =  PAM01325.0 
PVEL  =  IV(RHODS*B) 


— INITIALIZE  VARIABLES,  SET  FLAGS,  AND  CALCULATE  WET  BULB  TEMPERATURE _ 


CALL  TNITLAL 
CALLWETBULB 

- SPECIFIC  HEATS  AND  TRANSPORT  COEFFICIENTS  ARE  FOUND  FOR 

ESTIMATION  OF  DRYING  RATE  AND  CONTROL  VOLUME  LENGTH - 

CPA  =  CPAF<TAL/2.) 

CPF  =  CPFF(TPIV2.) 

HC  =  HCF(TAL,0.000 1  ,TPL) 

CALLRMODEL 
CPV  =  CPVFCITL/2.) 

YSAT  =  YWF(TWB) 

RKO  =  RKOF(TALYSAT.FRELTDEW,TWB) 

NV  =  RKO*(YSAT-YA) 

DELZ  =  (-L*DELTX)/(NV) 

- ESTIMATION  OF  DRYING  RATE  AND  CONTROL  VOLUME  LENGTH  FOR 

DIELECTRIC  DRYING - 

IF  (MWP.EQ.l)THEN 
PD  =  POWER/(B*A) 

NVO  =  POWER/(A*HFGF(100.)) 

DMF  =  1.+  HC*A*(TA-TSURF)/POWER 
QGEN  =  QGENF(PDEX)MF) 

NV  =  FREL*NVO*DMF 
DELZ  =  (-L*DELTX)/(NV) 

ENDIF 
C 

ENTHAL  =  CPAF(TAL/2.)*TAL+YAL*(HFGF(TDEW>+CPVF((TAL+TDEW)/2.)* 

*  (TAL-TDEW)+CPFF(TDEW/2.)*TDEW) 

ENTHO  =  ENTHAL 

C 

IF(ISO.EQ.l)THEN 

QAS  =  NV*(HFGF(TAL>+CPVF(TAL/2.)*TAL-CPVF(TPIV2,)*TPL>+ 

*  HC*(TAI^TPL) 

ENDIF 

C 

K  =  0 
C 

DO  100  I  =  l.N+l 
IF  (I.EQN+1)THEN 

XMC  =  XMCL+DELTX*FLOATG-l  )-DELTX/2. 


YA  =  YAL+DELTY*  rLOAT(I- 1  )-DELTY/2. 

DELTAX  =  DELTX/2. 

ELSEIF  (I.EQ.l)THEN 
XMC  =  XMCL 
YA  =  YAL 

DELTAX  =  DELTX/2. 

ELSE 

XMC  =  XMCL+DELTX*FLOAT(I- 1 ) 

YA  =  YAL+DELTY*FLOAT(I-l) 

DELTAX  =  DELTX 
ENDIF 

CALL  RMODEL 
CALL  ASSIGN  (X.O.) 

CALL  NLSYST  (FCNXFJJEQ.K) 

CALL  ASSIGN  (X,l.) 

ZPOS  =  ZPOS+DELZ 

NTU  =  NTU+DELTYANV/KKO) 

DRYTIM  =  ZPOS/PVEL 
IF  (YA.GT.YSAT)  IFLAG1  =  1 
IF  (TADSAT.GE.100.)  IFLAG2  =1 

TMEANO  =  (TSURF+TPV2. 

ENTHO  =  ENTHAL 
TPOLD  =  TP 
TSURFO  =  TSURF 

IF  (IP.EQJPRINT.ORJFLAG1  .EQ.l  OR.IFLAG2  .EQ.1)THEN 
IP=1 

YS  =  (YS  AT -YA)*FREL+YA 

- WRITE  VALUES  OF  VARIABLES  FOR  CURRENT  STEP 

CALL  WRITER  (1.) 

ELSE 
IP  =  IP+1 
ENDIF 

; - IF  AIR  IS  SATURATED,  EXIT  LOOP - 

IF  (IFLAG1  .EQ.l  .ORJFLAG2  EQ.l)  GO  TO  1 10 
100  CONTINUE 
1 10  CONTINUE 

: - WRITE  DRYER  LENGTH.DRYING  TIME  AND  NTU  — 

CALL  WRITER(2.) 

IF  (IFLAG1  EQ.l)  STOP  CONDENSATION  Y-AIR  >  Y-SURFACE' 
IF  (IFLAG2  .EQ.l)  STOP  ADIABATIC  TEMP.IS  G.E.100  C 
STOP  END  OF  PRGM.' 


£********************************************************************* 

c  - FCN  SUBROUTINE  CONTAINS  EQNS  SOLVED  BY  NLSYST 

£•••»•**•*•«*»••*•*•»•***•*••*••••••**•**••*•••*•****  ***************** 

SUBROUTINE  FCN  (X,F) 

DIMENSION  X(10),F(10) 

REAL  NVO,KDJLNV 
C 

COMMON/FUNCT/DELZ.DRYTIMENTHAUENTHOX»MFJ4EQjmj>rV(NVOJ>D, 

*  QAS,QCOND,QCONV,QEVAP,POWER.QGEN,QSENSD,QSENSW,TA,TADSAT, 

*  TDEW,TMEANO,TP,TPOLD,TSURF,TS  URFO.TWB.YS  AT .ZPOS.QV AP.QSOL 
COMMON  A3, D.DELTAX.GJJSO.UMWP.PWS.RELHUM.TAI.XMC.PA.YA 
common/propa:pa,cpai.cpf,cpfa,cps,cpv,cpva,fka,hfga,hfgs, 

*  HFGW,KD.RHODS,VISC 

COMMON/PMODEUEQUILF.SW,EXPON,FRACT,FREL,XCR,XNORM,SF,EFC,EFG. 

*  FI.XB.XCD.FF 
C 

C  -RENAME  SUBSCRIPTED  VARIABLES  WITH  THE  NAMES  USED  IN  THIS  CODE 

C 

CALL  ASSIGN  (X,l.) 

C 

IF  (ISO.EQ.l)THEN 
TAIR  =  TAI 
CPA  =  CPAI 
ELSE 

TAIR  =  X(2) 

CPA  =  CP  AF(T AIR/2.) 

ENDIF 

C 

CPF  =  CPFF(TP/2.) 

CPV  =  CPVF(TP/2.) 

CPVA  =  CPVF((TAIR+TDEWy2.) 

CPFA  =  CPFF(TDEW/2.) 

HFGW  =  HFGF(TP) 

HFGS  =  HFGF(TWB) 

HFGA  *  HFGF(TDEW) 

C 

C - COMPUTE  NORMALIZED  MOISTURE  CONTENT,  RELATIVE  DRYING  RATE, 

C  AND  THE  EVAPORATIVE  FRONT  POSITION - . 

C 

CALL  RMODEL 

HC  =  HCF(TAIR,NV,TSURF) 

C 

C - COMPUTE  THE  FLUXES  TO  BE  USED  IN  THE  ENERGY  BALANCE  EQUATIONS. 

C 

QVAP  =  NV*CPV*(TSURF-TP)/2. 

QSOL  =  NV*CP5*TP*EFG 
QCONV  =  HC*(TAIR-TSURF) 

QSENSW  =  L*((  1  .-FRACD*CPS+XMC*CPF)*(TP-TPOLD)/DELZ+QSOL 
QSENSD  =  FRACT*L*CPS*(((TSURF+TPy2.)-TMEANO)/DELZ-QSOL 
QEVAP  =  NV*HFGW 

RKO=  RKOF(TAIR,YSAT,FRELTDEW,TWB) 

YWSAT  =  YWF(TWB) 

IF  (FRACT.GT.O.)  QCOND  =  KD*(TSURF-TP)/(B*FRACr) 

C 

C - -COMPUTE  THE  HEAT  GENERATION  TERM. - 

C 


IF  (MWP.EQ.l)THEN 


nooo  n  non 
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IF  (FRACT  i.E.0.)  THEN 
DMF  =  1 .+  QCONV*  A/POWER 
ELSE 

DMF  =  l.+  QCOND*  A/POWER 
ENDIF 

QGEN  =  QGENF(PD.BX>MF) 

ELSE 
QGEN  =  0. 

ENDEF 

- SET  THE  GOVERNING  EQUATIONS  TO  BE  SOLVED.  F(l),  F(3).  F(4),  F(5), 

AND  F(6)  ARE  CONSTITUITIVE  EQUATIONS.  F(2X  F(7),  AND  F(8)  ARE  THE 
ENERGY  BALANCE  EQUATIONS. - 

F(l)  =  L*DELTAX  +  NV*DELZ 

F(2)  =  G*(ENTHAL-ENTHOyDELZ+QCONV<JAS-NV*(CPF*TP+HFGW) 

F(3)  =  CPA*TAIR+YA*(CPFA*TDEW+HFGA+CPVA*(TAIR-TDEW))-ENTHAL 
F(4)  =  YWF(TDEW>- Y  A 

F(5)  =  YWF(TADSAT)-(CPAF(TADSAT/1)*TADSAT-CPA'TAIR- 

*  YA*(HFGA+CPVA*TAIR-CPFF(TADSAT/2.)*TADSAT)y(CPFF 

*  (TADSAT/2.)*TADSAT-HFGF(TADSAT)-CPVF(TADSAT/2.)*TADSAT) 

IF  (MWP.EQ.O)  THEN 

F(6)  =  YWFfTWB)  -  YSAT 
ELSE 

F(6)  =  YWF(TADSAT)  -  YSAT 
ENDIF 

IF  (TPOLD.LT.  99 .9)  THEN 
IF  (FRACTDE.O.)  THEN 
F(7)  =  QCONV TJS  ENSW+QGENTJEVAP 
F(8)  =  TSURF-TP 
ELSE 

F(7)  =  QCOND-QSENSW+QGENTJEVAP 
F(8)  =  QCOND+QSENSD-QCONV  +  QVAP 
ENDIF 

IF  (MWP.EQ.O)  THEN 

F(9)  =  FREL*RKO*D*ALOG((D+YSATy(D+YA))-NV 

F(10)  =  RKO*ALOG((D+YWSATV(D+YA))*HFGS*D-HC*(TAIR-TWB) 

ELSE 

F(9)  =  FREL*NVO*DMF-  NV 
ENDIF 

- -IF  BOILING  POINT  DRYING  DEVELOPS,  THE  FOLLOWING  EQUATIONS  NEED  TO 

BE  USED. - 


ELSE 

IF  (XNORM  .LT.  1 .)  THEN 
LF(FRACT.GT.O.)THEN 

F(7)  =  QCONV -QSENSD+QGENjQEVAP  -  QVAP 
IF  (MWPEQ.l)  F(8)  =  FREL*NVO*DMF  -  NV 
IF  (MWP.EQ.O)  F(8)  =  QCOND/HFGW  -  NV 
ELSE 

IF  (MWP.EQ.l)  F(7>=  FREL*NVO*DMF  -  NV 
IF  (MWP.EQ.O)  F(7)=  QCONV/HFGW  -  NV 
ENDIF 

ELSEIF(FRACT.GT.O.)  THEN 
F(7)  =  QCONV43SENSD+QGEN-QEVAP  -  QVAP 


non  n  nnnn  n  n  nnnnnnnn 
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IF  (MWP.EQ.l)  F(8)=  FREL*NVO*DMF  -  NV 
IF  (MWP.EQ.O)  F(8>=  QCOND/HFGW  -  NV 
ENDIF 
ENDIF 
RETURN 
END 


-SUPPORTING  SUBROUTINES 

(********•************•************************************************ 


- FUNCTION  TO  COMPUTE  RELATIVE  DRYING  RATE  AND 

EVAPORATIVE  FRONT  POSITION - 

SUBROUTINE  RMODEL 

REAL  NVO.KD.LXEXENGTH.NUCl^UEXPl.NUL 

COMMON  A.B.D.DELT  AX, GJJSO.LMWP.PWS.RELHUM.TAI.XMC.PA.YA 
COMMON/PMODEL/EQUILF,SW,EXPON,FRACT.FREL,XCR,XNORM,SF,EFC,EFG, 

*  FDCB.XCD.FF 

COMMON/ITNrr/IFLAGlJFLAG2,K,YAL  XMCLXMCR.TPLTAUFLAG 
COMMON/FUNCT/DELZ,DRYTTM,ENTHALENTHOT)MFJ'IEQJJTUJIVJ^VOJ>D, 

*  qas.qcond,qconv,qevap,power,qgen,qsensd,qsensw,tatadsat, 

*  TDEW,TMEANO,TP,TPOLD,TSURF,TSURFO,TWB,YSAT.ZPOS,QVAP,QSOL 

IF  (EQUILF.GT.0.0)  then 
RELH  =  RHCT.YA) 

RELHUM  =  RELH 
XEQ  =  EQUILF*RELHUM 
IF  (XEQ.GT.XMC)  STOP  XEQ  >  XMC 
ELSE 
XEQ  =  0. 

ENDIF 


- COMPUTE  THE  NORMALIZED  MOISTURE  CONTENT  FOR  THE  FALLING  RATE 

AND  WARM-UP  PERIODS. - 

XNORM  =  (XMC-XEQy(XCR  -XEQ) 

XW  =  (XMC  -  XMCL)/(XB  -  XMCL) 

XF  =  (XMC-XMCR)/(XCD-XMCR) 

IF  (XNORM.LT.O.)  STOP  XNORM  IT.  ZERO' 

- THIS  BLOCK  DEALS  WITH  THE  FALLING  RATE  PERIOD - 

IF  (XNORM. LT.  1.0)  THEN 
IF  (MWP.EQ.l) THEN 
IF(SF.EQ.l)  FREL  =  FF  +  (1.-  FF)*XF 
IF(SF.EQ.2)  FREL  =  FF+2.*(1.-FF)*XF-(1.-FF)*XF**2 
IF(SF.EQ3)  FREL  =  FF+3.*(1.-FF)*XF-3.*(1.-FF)*XF**2+ 

*  (l.-FF)*XF**3 

FRACT  =  l.-(XMC/XMCL)**EFC 

EFG  =  (B*EFC/XMCL)*(XMC/XMCL)**(EFC-1) 

ELSE 

FREL  =  XNORM**EXPON 
FRACT  =  1  .-SQRT(XNORM/FREL) 

EFG  =  0. 


non 
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■THIS  BLOCK  DEALS  WITH  THE  WARM-UP  PERIOD - 

ELSEIF  (MWP.EQ.l)THEN 
IF(SW.EQ.l)  FREL  =  FI  +  (1.-  FI)*XW 
IF(SW.EQ.2)  FREL  =  FI  +  2.*(1.-  FI)*XW  -  (1.-  FI)*XW**2 
IF(SW.EQ.3)  FREL  =  FI+3.*(1.-FI)*XW-3.*(1.-FI)*XW**2+ 

*  (l.-FI)*XW**3 
IF  (FREL.GT.l.)  FREL=1. 

FRACT  =  l.-(XMC/XMCL)**EFC 
EFG  =  (-B*EFC/XMCL)*(XMC/XMCL)**(EFC-1) 

ELSE 
FREL=  1. 

FRACT  =  0. 

EFC  =  0. 

ENDIF 
RETURN 
END 

£**•**•***•••••*•*•**••••••****•*•••*****•*••••*•**•**•••••*••»•****** 

c - FUNCTION  TO  COMPUTE  THE  VOLUMETRIC  HEAT  GENERATION. - 

p************************ ************ ********************************** 

FUNCTION  QGENF(PD.B  J3MF) 

COMMON/PMODEIVEQUILF^W,EXPON,FRACr.FREUXCR^NORM,SF,EFC.EFG, 

*  FI,XB,XCD,FF 

IF  (XNORM.LT.  1)  THEN 
QGENF  =  PD*FREL*B/DMF 
ELSE 

QGENF  =  PD*B 
ENDIF 
RETURN 
END 
C 

£*************************************************************** 

C - -SUBROUTINE  TO  CALCULATE  INITIAL  GUESS  FOR  WET  BULB  AND 

DEW  POINT  TEMPERATURE - 

C*************** ********* *************************************** 

SUBROUTINE  WETBULB 
DIMENSION  F(16),X(16),XT(2),FT(2) 

REAL  KD.L.LEXENGTH  JWC1 JWEXP1  JVUL»NTUJTVJTVO 
EXTERNAL  TFCN 
C 

COMMON/FUNCT/DELZ.DR  YTIM^NTHALENTHOJDMFNrEQ^TU^VJTVOTD, 

*  QAS,QCOND,QCONV.QEVAP,POWER,QGEN,QSENSD,QSENSW,TA,TADSAT. 

*  TDEW,TMEANO,TP.TPOLD,TSURF,TSURFO,TWB,YSATZPOS,QVAP.QSOL 
COMMON  A3,D,DELTAX,G,USO,UMWP,PWS,RELHUM.TAI^MC,PA,YA 

C 

K  =  0 
NTQ  =  2 

YWS  =  YWF(TA) 

PWA  =  PA*YA/(DFYA) 

TA F  =  (TA*9V5.}+32. 

TWMIN  =  TAF*TAF*(-2.241468E-03)+.8315167*TAF-3.855549 
TGU  =  TWMIN-KTAF-TWMIN)*  PWA/PWS 
TGUC  =  «TGU-32.)*5./9.)*1.4 
XT(1)  =  TGUC 


onnnnn 
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XT(2)  =TGUC 
C 

CAUL  NLSYST  (TFCN.XT.FT,NTQ,K) 

C 

TWB  =  XT(2) 

TDEW  =  XT(1) 

TADSAT  =  XT(2) 

RETURN 

END 

**************************************************************** 

- SUBROUTINE  TO  RENAME  SUBSCRIPTED  VARIABLES  WITH  THE  NAMES 

USED  IN  THIS  CODE.  THE  VARIABLES  AND  NUMBER  OF  EQUATIONS 

MUST  MATCH  THOSE  OF  SUBROUTINE  FCN - 

**************************************************************** 

SUBROUTINE  ASSIGN  (X.T) 

DIMENSION  X(16) 

REAL  NVNVO 
C 

COMMON/FUNCT/DELZ,DRYTrM,ENTHAL,ENTHO.DMFJVEQ,NTUNV.NVOTD. 

*  QAS,QCOND,QCONV,QEVAP,POWER,QGEN,QSENSD,QSENSW,TATADSAT. 

*  TDEW.TMEANO.TP.TPOLD, TSURF, TSURFO, TWB, YSAT,ZPOS,QVAP,QSOL 
COMMON  A^,D,DELTAX,G,IJSO,UMWP,PWS,RELHUM,TAI^XMC,PA,YA 
COMMON/PMODEL/EQUI  LF,SW1EXPON,FRACT,FRELvXCR,XNORM,S  F.EFC.EFG, 

*  FLXB,XCD,FF 
C 

IF  (T.EQ.O.)  THEN 
IF  (TPOLD.LT.99.9)  THEN 
X(1)  =  NV 

IF  (ISOEQ.O)  X(2)  =  TA 
IF  (ISO.EQ.l )  X(2)  =  QAS 
X(3)  =  TP 
X(4)  =  DELZ 
X(5)  =  ENTHAL 
X(6)  =  TSURF 
X(7)  =  TDEW 
X(8)  =  TADSAT 
X(9)  =  YSAT 
IF(MWP.EQ.O)  THEN 
X(10)  =  TWB 
NEQ=  10 
ELSE 
NEQ  =  9 
ENDIF 
C 

ELSE 

IF  (ISO.EQ.O)X(l)  =  TA 
IF  (ISO.EQ.l)  X(l)  =  QAS 
X(2)  =  DELZ 
X(3)  =  ENTHAL 
X(4)  =  TDEW 
X(5)  =  TADSAT 
X(6)  =  YSAT 
IF  (XNORM.LT.l.)THEN 
IF  (FRACT.GT.0.)  THEN 
X(7)=  TSURF 
X(8)  =  NV 
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X(7)  =  NV 
NEQ=7 
ENDIF 

ELSEIF(FRACT.GT.O.)THEN 
X(7)  =  TSURF 
X(8)  =  NV 
NEQ  =  8 
ELSF 
NEQ  =  6 
ENDIF 
ENDIF 

ELSEIF  (TPOLD.LT.99.90)  THEN 
NV=X(1) 

IF  (ISO.EQ.O)  TA=X(2) 

IF  OSO.EQ.l)  QAS=X(2) 
TP=X(3) 

DELZ=X(4) 

ENTHAL=X(5) 

TSURF=X(6) 

TDEW=X(7) 

TADSAT=X(8) 

YSAT=X(9) 

IF  (MWP.EQ.O)  THEN 
TWB  =  X(10) 

ENDIF 

ELSE 

IF  (ISO.EQ.O) TA  =  X(l) 

IF  (ISOfQ.l)QAS  =  X(l) 
DELZ  =  X(2) 

ENTHAL  =  X(3) 

TDEW  =  X(4) 

TADSAT  =  X(5) 

YSAT  =  X(6) 

IF  (XNORM.LT.l.)  THEN 
IF  (FRACT.GT.O.)  THEN 
TSURF  =  XC7) 

NV  =  X(8) 

ELSE 
NV  =  X(7) 

ENDIF 

ELSEIF(FRACT.GT.O.)  THEN 
TSURF  =  X(7) 

NV  =  X(8) 

ENDIF 

ENDIF 

RETURN 

END 


uuu  uuu 
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£«*********•***••*****•*********•*******•***********»****•********• 

SUBROUTINE  WRITER  (PEN) 

REALNTU.NV.NVO 
DIMENSION  X(16) 

C 

COMMON/FUNCT/DELZ.DRYTTM.ENTHAL^NTHO.DMF.NEQ.NTU.NV.NVO.PD, 

*  QAS,QCOND,QCONV,QE V AP.POWER.QG EN.QSENS D,QS ENS W.TA.TADS AT, 

*  TDEW,TMEANO,TP,TPOLD,TS U  RF.TS  URFO,TWB, YS AT ,ZPOS ,Q V AP.QSOL 
COMMON  A3.D,DELTAX,GJJS0,LMWP,PWS,RELHUM,TAI^CMC,PA.YA 
common/pmodeivequtlf,sw,expon,fract,frel>xcr.xnorm,sf,efc,efg. 

*  FLXB.XCD,FF 

COMMON/CONV/ATM,BETACPRAT,CPY,DAB.HC.LE,LENGTILNUCl.NUEXPl.NUL. 

*  PHIH,PR,PSYRAT,REURHOA,RHOMA.RKO,SC,VEL,YDAYDS,YS 

- -WRITE  HEADERS  FOR  OUTPUTFILES . - 

IF  (PEN.EQ.O.)  THEN 
WRITE  (1.10) 

10  FORMAT(2X,TIME',3X,'ZPOS',5X,X3X,‘YA/YSAT',3X,'DMF,TSURF, 

*  4X.TF.5X,-NV,5X;FREL') 

WRITE  (2,20) 

20  FORMAT(3X,T,5X.'XMC',9X,TA',9X,’YA'.9X,TP‘,9X,'YSAT,,9X.,RH') 

WRITE  (3.85) 

35  FORMAT  (1X,'XMC,8X,'NV',7X.'QGEN\8X, 

*  T)MF,6X,'FREL') 

. . WRITE  VARIABLES  FOR  EACH  STEP . 

ELSEIF  (PEN.EQ.l.)THEN 

WRITE(1,120)  (DRYTTM/60.),ZPOS,XMC,(YA/YSAT),DMF,TSURF, 

*  TP.NV.FREL 

120  FORMAT(1X,3(F7.2,'.'),F4.2,',',3(F6.1,’,’),E10.3,',',F7.3) 

WRITE  (2,70)  LXMC.TA.YA.TP.YSAT.RELHUM 
70  FORMAT  (1X.I3.7,F7.4,7,F7.4,7,F8.4,7.F8.4,7,F7.4. 

*  ,,',F7.4,',,tF6.1,7.F6.1.',',E10.3,',’.F5.2,’,',F5.2) 

WRITE  (3.95)  XMC,NV,QGEN,DMF,FREUNVO 
95  FORMAT  (1X.9(E12.6,',')) 

WRITE  (430)  (DRYTIM/60.)XMC.TA,TP,TADSAT.TSURF,TAI 
30  FORMAT  (1X,7(7F93)) 

WRITE  (535)  XNORM,XMC,(NV*7920.*A) 

35  FORMAT  (1X3(F73,7),E93,7) 

WRITE  (7.95)  ZPOS.HC.RKO.PS  YRAT 

WRITE(8,80)  XMC.QCOND.QSENSW.QEVAP.QSENSD.QCONV, 

*  QAS.QGEN.QVAP.QSOL 
80  FORMAT  (1X,10(F63,',,)) 

WRITE  (10,40)  ZPOS,YA,YSAT,(YSAT-YA) 

40  FORMAT  (FI  1.4,7,F10.8,7,F10.8,7,F10.8) 

WRITE  (1 1 30)  LXMC,YA,TWB,YSAT,(YSAT-YA)JW 
50  FORMAT  (I3,5(7,F11.7),E10.5) 

WRITE  (12,60)  TA.YATDEW.TADS  AT.TWB.QAS 
60  FORMAT  (1X,6(F6.2,7)) 

WRITE  (13.*)  PR,SC,LE,PHIH 

WRITE  (14,*)  BETA.CP  RAT.LE.PS  YRAT 

WRITE  (15,*)  YAENTHAL,  FRA  CT 

WRITE  (16,*)  ZPOS.TAQCONV.QEVAP.QSENSW.QAS 


non  non  n o n  noon 


WRITE  (17,90)  XMC,YAJCNORM,FREL,NV,TSURF,ZPOS 
90  FORMAT  (F7.4,7.F7.5.7.F7.5,7,F5.3.7,E12.4,7,F7.2, 

•  '.'.F10.3) 

- WRITE  THE  DRYER  LENGTH,  DRYING  TIME,  NTU,  DIELECTRIC  FLUX, 

AND  LATEST  AIR  TEMPERATURE. - 


ELSE 

WRITE  (1,77)  ZPOS,(DRYTIM/60.),ABS(NTU),(QGEN*DMF/FREL),TA 
77  FORMAT  (/.1X/DRYER  LENGTH  =  ’,Fl  1.2.’  METERS’,4X,  DRYING  TIME 

*  -'.F12.Z'  MINUTES’.//  NUMBER  OF  TRANSFER  UNITS,  NTU  =’. 

*  F7.2,4X,’GENE RATED  POWER  =’,F7.2.’  KW/SQ.METER7.  AIR  TEMP.=', 

*  F7.1,'  DEGREES  C) 

ENDIF 

RETURN 

END* 

***********************++**♦♦**+**#*++*♦+♦♦♦♦*+***+#+*****+***+*•#*»* 

- FUNCTION  TO  COMPUTE  SPECIFIC  HEAT  OF  AIR - 

***************************************************  ***************** 

FUNCTION  CPAF  (D 
TK  =  273.15+T 
TKl  =  TK/100. 

CPAF1  =  (.2l9+.342*(TK*9./5.)/lE4-.293*(TK*9./S.)**2/lE8)*4.1869 
CPN2  =  39.06-5 1 2.79/(TKl  *SQRT(TK  1  ))+1072.7/TK  1  **2-820.4/TKl  **3 
CP02  =  37.432+0.020102 *TK1*SQRT(TK1)-178.57/(TK1*SQRT(TK1))+ 236. 88 

*  /TK1**2 
CPAR  =  .5203 

CPAF2  =  (.78 1  *CPN2+.2 1  *CP02+.009*CPAR)/28 .9645 
CPAF  =  CPAF2 
RETURN 
END 

******** •*•**«•*«»•*******•**••****•*****» ************************** 

- FUNCTION  TO  COMPUTE  SPECIFIC  HEAT  OF  SAT  LIQUID  WATER 

******************************************************************** 

FUNCTION  CPFF  (T) 

C20  =  1.02493 
C2I  =  (-4.149E-04)*2. 

C22  =  (3.07768E-06)*3. 

C23  =  (-1.2606E-08)*4. 

C24  =  (3.0658 IE- 11  )*5. 

C25  =  (-3.843E-14)*6. 

C26  =  (1.9907E-17)*7. 

CPFF  =  (C20+C21*(T*9./5.+32.)+C22*(T*9./5.+32.)**2+C23*(T*9./5. 

*  +32.)**3+C24*(T*9./5.+32.)**4+C25*(T*9./5.+32.)**5+C26*(T*9. 

*  /5.+32.)**6)*4.1869 
RETURN 

END 

****************************************************************** 

- FUNCTION  TO  COMPUTE  SPECIFIC  HEAT  OF  WATER  VAPOR  - 

****••**••*•••»••**•••**********%***** **************** ************* 
FUNCTION  CPVF(T) 

TK  =  273.15+T 
TKl  =TK/100. 

IF  (TK.LT.0.0)  STOP  TEMP  <  0.0  DEG  CPVF 

CPVF  =  (143 .05-1 83. 54*TK1**0. 25+82, 75 1*SQRT(TK1)-3.6989*TK1)/ 

*  18.01534 


-Ci "o'.v.v.v; 
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RETURN 

END 

^***»*«**»»*»**«****»****. *•••••»************•*******•*******••*•*** 

C  -  FUNCTION  TO  COMPUTE  SATURATED  HUMIDITY  - 

C************* . . . . . ***...*.********.. 

FUNCTION  YWF  <T) 

REAL  NVO.KD.ULEXENGTRNUC  1  NUEXP  1  .NUL 
C 

COMMON  A,B.D,DELTAX,GJ.ISO.L.MWP,PWS,RELHUM,TAI.XMC,PA,YA 
C 

TK  =  273.15+T 
TK  =  ABS(TK) 

IF  (T.LT.0.0)  STOP  EXPLODING  SOLUTION.  T  <  0  K  YWF 
C8  =  -5800.2206 
C9  =  1.3914993 
CIO  =  -0.04860239 
Cll  =  0.4176478E-04 
C12  =  -0.14452093E-07 
C13  =  6.5459673 
C 

PWS  =  EXP(C8/TK+C9+C10*TK+C11*TK**2+C12*(TK)**3+C13*ALOG(TK)) 

IF  (PA.LE.PWS)  THEN 
YWF  =  9999999.9 
ELSE 

YWF  =  ,62198*PWS/(PA-PWS) 

ENDIF 

C 

RETURN 

END 

£**•****•*•*•••**••*•*•••******•**••*••*•*•*•*••*••**•*•***•*••*••**•• 

c - FUNCTION  TO  CALCULATE  LATENT  HEAT  OF  VAPORIZATION 

£*******»**•*»*««*»** ****** ******************************************** 

FUNCTION  HFGF  (D 
TAF  =  T*9./5.+32. 

TAF  =  Auov  ‘  AF) 

IF  (TAF.LT.0.0)  STOP  EXPLODING  SOLN,  TA  <  0.  HFGF 
HFGF  =  (1094.40516*EXP( -5.3882271 1E-04*TAF))*2.326 
RETURN 
END 

£*«*********************«********************************»*******»**** 

c - FUNCTION  TO  CALCULATE  RELATIVE  HUMIDITY 

£•*****•**••*••*•*•*•*•*•*••**»*»**••*»•••**•••****•****•*•*•****•**•* 

FUNCTION  RH  (TN.YA) 

REAL  NVO.MU 
IF  (TN.GT.99.5)  THEN 
RH  =  YA/(YA+1.) 

ELSE 

YSAT  =  YWF(TN) 

MU  =  YA/YSAT 

PRR  =  YSAT/(.62198+YSAT) 

RH  =  MU/(1.-(1.-MU)*PRR) 

ENDIF 

RETURN 

END 
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c**»* . ************ . 

FUNCTION  DHFGF  (TN) 

T  =  TN*9./5.+32. 

DHFGF  =  (-.59690355*EXP(-5.38822?11E-04*T))*4.1868 
RETURN 
END 


- FUNCTION  TO  COMPUTE  CONVECTIVE  HEAT  TRANSFER  COEFF.  - 


FUNCTION  HCF  (T,NVS,TS) 

REAL  NVOJCD.ULE^LENGTHNUCl^fUEXPl^IUUNVS 
C0MM0N/PR0P/CPA.CPAI.CPF,CPFA,CPS,CPV,CPVA,FKA.HFGA.HFGS, 

*  HFGW,KD,RHODS,VISC 

COMMON/CONV/ATM^ETACPRAT,CPY.DAB,HC^ELENGTHNUCl,NUEXPl.NUL 

*  PHIH,PR,PSYRAT,RELRHOA,RHOMA,RKO,SC,VEL.YDAYDS,YS 

C 

TF  =  (T+TS)/2. 

VISC  =  VISCFCTF) 

REL  =  VEL* LENGTH/VI  SC 
NUL  =  NUC1*ABS(REL)**NUEXP1 
HC  =  NUL*  FKAER(TF)/LENGTH 
E  =  NVS*CPVF((TF)/1)/HC 
IF  (E.GT.1.E-9)THEN 
PHIH  =  E/(EXP(E)-1.) 

ELSE 

PHIH  =  1.00 
ENDIF 

HCF  =  HC*PHIH 
HC  =  HCF 
RETURN 
END 

q** **********************************************  ******************* 

C  —  FUNCTION  TO  COMPUTE  THERMAL  CONDUCTIVITY  OF  AIR - 

£***********«*****»**********•***•»»»***»»*»*»*»•*»*»*»**»**»»*»»**« 

FUNCTION  FKAIR  (T) 

FKAIR  =  (24.13022+  ,0749196*T-3.635783E-05*T**2+2.32043E-08*T** 

*  3)*l.E-06 
RETURN 
END 

£**•***•**••****•••»*****•••*•* ****************************************** 

C -  FUNCTION  TO  COMPUTE  AIR  VISCOSITY - 

C************************************************************************ 

FUNCTION  VISCF  (T) 

VISCF  =  (13. 28319+8. 803  2E-02*T+9.883066E-05*T**2-3.1548E-08*T**3)* 

*  l.E-06 
RETURN 
END 

£**************************************»********»******** ****** ******** 

C - FUNCTION  TO  COMPUTE  MASS  TRANSFER  COEFFICIENT - 

Q*** ***********************************************************  ******** 

FUNCTION  RKOF  (TA,YW,FC,TDP,TWB) 

REAL  KO.MW.MA.MBS.MBA 
REAL  NVO,KD.LJLELENGTH,NUC  1 J9UEXP 1 .NUL 
C 

COMMON  AB,D,DELTAX,GJ,ISO,LAIWP,PWS1RELHUM,TAIJ(MC,PAYA 
COMMON/PROP/CPA, CPAI,CPF,CPFA,CPS,CPV,CPVA,FKA,HFGA,HFGS, 

*  HFGW,KD,RHODS,VISC 


COMMONATONV/ATM.BETA,CPRAT,CPY,DAB,HC.LE,LENGTH.NUCl,NUEXPl.NUL, 
*  PHIH.PR.PSYRAT.RELRHOA.RHOMA.RKO.SC.VEUYDA.YDS.YS 

YS  =  (YW-YA)*FC+YA 
TF  =  (TA+TWBV2. 

YF  =  (YA+YWV2. 

IF  (Y  A.LT.  YW)  THEN 
MW  =  .01801534 
MA  =  .0289645 
YDS  =  YW/(D+YW) 

YDA  =  YAAD+YA) 

MBS  =  MW*YDS+MA*(1.-YDS) 

MBA  =  MW*YDA+MA*(1.-YDA) 

BETA  =  MBS*  ALOG(MB  S/MB  A  )/((MW  ■  MA)*D*  ALOG((D+ YWy(D+Y  A))/(D+YW)) 
ELSE 

BETA  =  1.01 
ENDIF 


TK  =  273.15 

DAB  =  1.1998806E-09*SQRT(TF+TK)*(TF+TK)*(TF+TK)**(.25) 

TCA  =  647.3 
TCB  =  132.48 
PCA  =  218.01 
PCB  =  37.10 
RMA  =  18.01534 
RMB  =  28.9645 

DENOM  =  (PCA*PCB)**(1./3.)*(TCA*TCB)**(5./12.)*SQRT(1./RMA+17RMB) 
RAD  =  SQRT(TCA*TCB) 

Cl  =  (l./RAD)**2.334 
C  =  3.64E-»*Cl*DENOM 
DAB1  =  C*(TF+TK)**2.334*l.E-4 

RHOA  =  353.13/(TF+TK) 

ATM  =  PA/101325. 

YWDP  =  YA 

DPP  =  YWDP/(D+YWDP) 

RHOMA  =  RHOA*(ATM-.3783*DPP) 

FKA  =  FKAIR(TF) 

VISC  =  VISCF(TF) 

CPA  =  CPAF(TF) 

CPY  =  CPA+YF*(CPVF(TFHDHFGF(TDP)) 

CPRAT  =  CPY/(1.+YF) 

PR  =  CPA*RHOA*  VISC/FKA 
SC  =  V1SC/DAB 
LE=  PR/SC 

PSYRAT  =  BETA*LE**(273.yCPA 
KO  =  HC*PSYRAT 
RKOF  =  KO 


RETURN 


9  7 


£•*****•••***••**••••••*•**••********•**•**•***•«**•*•****»•**•»****• 

C - SUBROUTINE  CONTAINING  EQUATIONS  FOR  DEW  POINT  AND  WETBULB 

C  TEMPERATURES.  THIS  SET  OF  EQS  IS  USED  INTHE  INITIAL  CALL  TO  NLSYST. 

£**•******•••******•**•****••**•*•*»*•*«*•*********»•**•*•*•***•*••** 

SUBROUTINE  TFCN  (X.F) 

DIMENSION  F(2),X(2) 

REAL  NVO.KD.ULE,  LENGTH  .NUC 1  NUEXP1.NUL 
C 

COMMON  A,B,D,DELTAX,GJ.ISO,UMWP,PWS,RELHUM.TAI,XMC.PA,YA 
C 

PWA  =  PA*YA/(D+YA) 

YWDP  =  YWF(X(1)) 

PWSDP  =  PWS 
YWSAT  =  YWF(X(2)) 

PWSAT  =  PWS 

HC  =  HCF(TAI,0.0,X(2)) 

RKO  =  RKOF(TAI,YWSAT,1.0,X(1),X(2)) 

HFGS  =  HFGF(X(2)) 

C 

F(l)  =  PWA-PWSDP 

F(2)  =  RKO*ALOG((PA-PWA)/(PA-PWSAT))*HFGS*D-HC*(TAI-X(2)) 

C 

RETURN 

END 

£*•***•*****•**••***•**•********••****•••*•****•*****•**,•*,»***•,•, 

c - -SUBROUTINE  THAT  SETS  FLAGS  AND  INITIALIZES  VARIABLES. — 

C****************************************** ****** ******************* 

SUBROUTINE  TNITIAL 

COMMON/FUNCT/DELZ,DRYTIM.ENTHAIJENTHOX>MFNEQ.NTUNV.NVO.PD, 

*  QAS.QCOND,QCONV.QEVAP,POWER,QGEN,QSENSD,QSENSW.TA,TADSAT, 

*  TDEW,TMEANO,TP.TPOLD,TSURF,TSURFO,TWB,YSATZPOS,QVAP,QSOL 
COMMON  A.B.D.DELTAX,GJJS0,UMWP,PWS,RELHUM,TAI.XMC,P  A.  YA 
COMMON/nNIT/IFLAGlTFLAGZK.YAUXMCL.XMCR.TPL.TAUIFLAG 
COMMON/PROP/CPA,CPAI,CPF,CPFA,CPS,CPV,CPVA,FKA,HFGA,HFGS. 

*  HFGW.KDJIHODS.VISC 
COMMON/NLSYS/MAXIT, DELTA  JCTOLFTOL 
REALNTU 

IFLAG  =  0 
IFLAG1  =0 
IFLAG2  =  0 
MAXIT=120 
XTOL  =  0.00001 
FTOL  =  0.0003 
K  =  0 

ZPOS  =  0.0 
NTU  =  0.0 
YA  =  YAL 
XMC  =  XMCL 
TSURF  =  TPL 
TA=TAL 
TAI-TAI 

CPAI  =  CPAF(TAL/2.) 

TMEANO  =TPL 
TP  =  TPL 
TPOLD  =  TPL 
TSURFO  =TPL 
RETURN 
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£**»**•**•*•****•***••**•**•******•*••***••»•*•*•*•**»•*****•****»•**»•*• 

C - NON-LINEAR  EQUATIONS  SOLVER.  TAKEN  FROM  REFERENCE  [15].-— 

£•****#****•*********•****•«****************»**»***********«*****•**«**** 

SUBROUTINE  NLSYST  (FCN.X.FNJ) 

DIMENSION  X(N),F(NXA(  10. 1 1 ).XSA  VE(1 1  ).FS  AVE(  1 1 ) 
TOMMON/NLSYS/MAXITDELTA.XTOL.FTOL 
C 

LOGICAL  PRINT 
C 

NEQ1  =  N 

IF  (N.LT.2.0RN.GT.15)  GO  TO  180 
PRINT  =  TRUE. 

IF  (I.NE.0)  PRINT  =  .FALSE. 

C 

NP  =  N+l 

Doi2orr=iMAxrr 
DO10IVBL=lN 
XSAVE(IVBL)  =  XflVBL) 

10  CONTINUE 
C 

CALL  FCN  (X.F) 

C 

ITEST  =  0 
DO  20 IFCN  =  IN 

IF  (ABS(F(IFCN)).GT.FTOL)  ITEST  =  ITEST+1 
FSAVE(IFCN)  =  F(IFCN) 

20  CONTINUE 
C 

IF  (NOT. PRINT)  GO  TO  50 
WRITE  (6,30)  nVX 

30  FORMAT  (IX.' AFTER  ITERATION  NUMBER'J3,'  X  AND  F  VALUES  ARE  'J 
*  .1X.13F10.4) 

WRITE  (6.40)  F 
40  FORMAT  (1X.13F10.4) 

C 

50  IF  (TTEST.NE.0)  GO  TO  60 
1  =  2 

RETURN 

C 

60  DO  80  JCOL  =  IN 

X(JCOL)  =  XSAVE(JCOL>+DELTA*XSAVE(JCOL) 

CALL  FCN  (X.F) 

DO70IROW=lN 

A(IROW,JCOL)  =  (F(IROW)-FSAVE(IROW))/(DELTA*XSAVE(JCOL)) 

70  CONTINUE 
C 

X(KX)L)  =  XSA  VE(JCOL) 

80  CONTINUE 
C 

DO90IROW  =  IN 
AOROWNP)  =  -FSAVEdROW) 

90  CONTINUE 

CALLEUM  (A.NNP.10) 

C 

DO  100 IROW  =  IN 


i 
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IF  (ABS(A(IROW,IROW)).LE.  1  .E- 14)  GO  TO  160 
100  CONTINUE 
C 

1TEST  =  0 
DO  110IVBL=  1.N 
X(TVBL)  =  XSAVE(IVBL)+A(TVB  UNP) 

IF  (ABS(A(IVBUNP)).GTJ<TOL)  ITEST  =  ITEST+1 
110  CONTINUE 
C 

IF  (ITEST  £Q.0)  GO  TO  130 
120  CONTINUE 
I  =  -l 
RETURN 
C 

1301=  1 

IF  ( JMOT.PRIND  GO  TO  1 50 
WRITE  (6,140)  IT,X 

140  FORMAT  (IX,' AFTER  ITERATION  NUMBER'43, 

*  'X  VALUES  (MEETING  XTOL)  ARE  V.1X.13F10.4) 

150  RETURN 
C 

1601  = -2 
WRITE  (6,170) 

170  FORMAT  (IX, 'CANNOT  SOLVE  SYSTEM.  MATRIX  NEARLY  SINGULAR') 

IF  (I.EQ.-2)  STOP  "NEARLY  SINGULAR  MATRIX1 
RETURN 
C 

1801  = -3 
WRITE  (6,190)  N 

190  FORMAT  (IX, "NUMBER  OF  EQUATIONS  PASSED  TO  NLSYST  IS  INVALID.  MU 
*ST  BE  1  <  N  <  13.  VALUE  WAS  J3) 

RETURN 

END 

£***•***» *♦********♦******♦*****•******+******•**************** 

SUBROUTINE  ELIM  (ABJ4J4PJ4DIM) 

DIMENSION  AB(NDIMNP) 

C 

NM1  =N-1 
DO  60 1  =  1JNM1 
C 

IPVT  =  1 
IP1  =  1+1 
DO  10  J  =  EP1.N 

IF  (ABS(AB(IPVTJ)).LTABS(AB(J.I)))  IPVT  =  J 
10  CONTINUE 
C 

IF  (IPVT  £Q  J)  GO  TO  30 
DO  20  JCOL  =  I.NP 
SAVE  =  AB(IJCOL) 

AB0JCOL)  =  AB  (IPVT  JCOL) 

AB  (IPVT  JCOL)  =  SAVE 
20  CONTINUE 
C 

30  DO  50  JROW  =  IP1.N 

IF  (AB(JROWJ).EQ.O)  GO  TO  50 
RATIO  =  AB(JROWJ)/AB(IJ) 

DO  40  KCOL=  IP1J4P 


AB(JROW.KCOL)  =  AB(mOW,KCOL)-RATIO*AB(l.KCOL) 

40  CONTINUE 
50  CONTINUE 
60  CONTINUE 

IF  (ABS(AB(NJ4)).LT.l  .E-14)  GO  TO  100 
NP1  =  N+l 

DO90  KCOL  =  NP1.NP 
AB(N,KCOL)  =  AB(N,KCOL)/AB(N,N) 

DO  80  J  =  2.N 
NVBL  =  NP1-J 
L  =  NVBL+1 

VALUE  =  AB(NVBL,KCOL) 

DO  70  K  =  UN 

VALUE  =  VALUE-AB(NVBUK)*AB(KJCCOL) 

70  CONTINUE 

AB(NVBUKCOL)  =  VALUE/AB(NVBUNVBL) 

80  CONTINUE 
90  CONTINUE 
RETURN 

100  WRITE  (6.110) 

110  FORMAT  (IX. 

*  SOLUTION  NOT  FEASIBLE.  A  NEAR  ZERO  PIVOT  WAS  ENCOUNTERED  •) 
WRITE  (6,*)  N,(AB(N,N)) 

IF  (ABS(AB(N.N)).LT.  1  .E- 14)  STOP  NEAR  ZERO  PIVOT 
RETURN 


no  n  n  on  n  r>  r>  n  n  n  n  n  n  n  n  n  n  n 

!  *  «  «- 


PROGRAM  NORM(INPUT, OUTPUT, DAT ATAPE1=DATA,TAPE9=  INPUT, 

*  T APE2.T  APE3  , TAPE4, TAPES  ,T APE6.T APE7.T APE8.TAPE  1 0.TAPE  1 1 ) 

*********************************************************************** 


NORMALIZING  CODE 

THIS  CODE  REQUIRES  NV  VERSUS  X  AND  TP  VERSUS  X  DATA  AS  INPUT. 

THE  OUTPUT  FILE,  WHEN  PLOTTED,  GIVES  THE  CHARACTERISTIC  DRYNG 
CURVE  FOR  THE  PRODUCT. 

THE  ACCOMPANYING  VARIABLE  DICTIONARY  PROVIDES  INFORMATION 
ON  VARIABLES  AND  SUBROUTINES  FOR  USE  OF  THE  CODE. 

i,**  *******************  •  **********  ***************************  *********** 


COMMON/CO  NV/ATM.B  CPRAT,CPY,DAB,HC,LE,LENGTH,NUC1,NUEXP1,NUL> 

*  PHIH.PR.PS  YRAT,REL,RHOA,RHOMA,RKO,SC,VEL,  YD  A,  YDS.YS 

DIMENSION  XMC(150),NV(150),TP(150),TSURF(150),NVD(150),XSTAR(150) 

REAL  NVO,KD.L,LE,LENGTILNUCl  JJUEXP1  JWUNVJ<VD 

READ  (9,*)  A3.POWER  JI.EFC JCB ,TAIR 

READ  (9,*)  VELJNUC 1 JWEXP1  ,L£NGTH,KD JPRINT 

- READ  DATA  OBTAINED  FROM  A  BENCH  SCALE  DRYING  TEST  ON  A  SMALL 

QUANTITY  OF  THE  PRODUCT.  THIS  CODE  ASSUMES  DRYING  RATE  IS  IN 
GRAMS  PER  MINUTE - 

READ  (1,*)  (XMCG).I=1J'I),(NVD(I),I-1,N),(TSURF(I)J=1,N),(TP(I),I=1J9) 


IP  =  IPRINT 

- NORMALIZE  EACH  DATA  POINT - 

DO  10  I=1,N 

NV(I)  =  NVD(iy60000. 

XSTAR(I)  =  (XMC(I)  -  XMC(1  ))/(XB  -  XMC(1 )) 

NVO  =  P/HFGF(TP1(I)) 

HC  =  HCF(TP1(I)JTV(D,TSURF(I)) 

FRACT  =  1 (XMC(I)/XMC(  I  ))**EFC 
QCONV  =  HC*A*(TAIR  -  TSURF(D) 

QCOND  =  *A*KD*(TSURF(I)  -  TP(I)V(FRACT*B) 

C 

C - DETERMINE  IF  AN  EVAPORATIVE  FRONT  IS  PRESENT  TO  CALCULATE 

THE  DRYING  MODE  FACTOR - 

C 

IF  (FRACT  .LT.0.02)  THEN 
DMF  =  1 .  +  QCONV/POWER 
ELSE 

DMF  =  1 .  +  QCOND/POWER 
ENDIF 
C 


c - COMPUTE  THE  RELATIVE  DRYING  RATE  AND  WRITE  THE  VALUES - 

C 

FREL  =  NV(R/(NVO*DMF) 

IF  (IP.EQJPRINT)  THEN 
IP  =  1 

WRITE(  2, 50)  XMC(I),FREL,FRACT,DMF 
WRITE  (3,50)  XMC(I),QCONV,QCOND,TSURF(I),FRACT 
WRITE  (5,50)  XMC(I),XSTAR(I),NVD(I),NV(I),FREL 
WRITE  (6,50)  XMC(I),NVD(I).TP(I) 

50  FORMAT(1X,7(F10.4,'.')) 

ELSE 

IP  =  IP+1 
ENDIF 

10  CONTINUE 
STOP 
END 

£******•*•*••••***« ****************************** *•**»****«>  *********** 

C - -SUPORTING  FUNCTIONS.  SAME  AS  IN  DRYER  CODE 

*********************************************************************** 

C 

FUNCTION  HFGF(D 
TAF  =  T*9./5.+32. 

TAF  =  ABS(TAF) 

IF  (TAF.LT.0.0)  STOP  -EXPLODING  SOLN,  TA  <  0.HFGF 
HFGF  =  (1094.405 1 6*EXP(-5. 388227 1 1  E-04*TAF))*2.326 
RETURN 
END 

£*••**••••••••••*•••*••••••**•*•**••**••••*••••*••*••••**•**••••••**•• 

FUNCTION  HCF  (TJ9VS.TS) 

REAL  NVO.KDUUEJ-ENGTUNUC  1 JNUEXPl  .NUUNVS 

COMMON/CONV/ATM3ETA.CPRAT,CPY,DAB,HC.LE,LENGTHJ4UClJ9UEXPU9UU 

*  PHIH.PR,PSYRAT.RELRHOA,RHOMA,RKO,SC,VEL,YDA.YDS.YS 
C 

TF  =  (T+TS)/2. 

VISC  =  VISCF(TF) 

REL  =  VEL*LENGTH/VISC 
NUL  =  NUC1*ABS(REL)**NUEXP1 
HC  =  NUL*FKAIR(TFyLENGTH 
E  =  NVS*CPVF((TF)/2.VHC 
IF  (E.GT.  1  .E-9)  THEN 
PHIH  =  E/(EXP(E)-1 .) 

ELSE 

PHIH  =  1.00 
ENDIF 

HCF  =  HC'PHIH 
C 

RETURN 

END 

c************************************************************************* 

FUNCTION  VISCFCT) 

VISCF  =  (13.283 19+8.8032E-02*T+9.883066E-05*T**2-3.1548E-08*T*  *3)* 

*  l.E-06 
RETURN 
END 


c 


FUNCTION  FKAIR  (T) 

FKAIR  =  (24.13022+.0749196*T-3.635783E-05*T**2+2.32043E-08*T** 

*  3)*l.E-06 
RETURN 
END 

£•***•******•*•••**•***•*•********»»•*••*••••*•»•••»••****••****•**• 

FUNCTION  CPVF(T) 

TK  =  273.15+T 
TK1  =TK/100. 

IF  (TK.LT.0.0)  STOP  TEMP  <  0.0  DEG  CPVF 

CPVF  =  (143.05-183.54*TK1**0.25  +82.751*SQRT(TK1)-3.6989*TK1)/ 

*  18.01534 
RETURN 


* 


B.l  HF.AT  AND  MASS  BALANCES 


The  control  volume  around  the  dielectric  dryer  was  divided  into  three 
control  volumes  to  represent  the  distincs  regions  inside  the  dryer  and  the  product. 
The  three  control  volumes  correspond  to  the  air  region,  the  dry  product  region,  and 
the  wet  product  region.  The  following  assumptions  were  used  to  perform  the 
analysis: 

•  The  process  is  steady-state,  and  the  following  conservation  laws  apply  where  M 
and  E  represent  mass  and  energy  respectively: 

Min  +  MSource  =  Mout  and  Ejn  +  Esource  =  Eout 

•  No  mass  source  is  considered  in  this  analysis. 

•  A  uniform  temperature  profile  exits  in  the  wet  region  represented  by  Tp. 

•  A  lumped  parameter  analysis  may  be  used  to  represent  thermal  and  mass 
transport  in  the  y-direction. 

•  Dryer  walls,  and  the  top  and  bottom  dryer  boudaries  are  insulated  and 
impermeable. 

•  Formulation  was  done  on  a  per  unit  of  dryer  width  basis. 

•  Dry  product  is  transparent  to  electromagnetic  energy. 

•  No  heat  conduction  in  the  z  and  x  direction 

•  No  product  is  entrained  in  the  air  stream. 

•  The  solid  and  air  mass  flow  rates  are  considered  constant 

•  A  linear  profile  exists  between  the  surface  temperature, Ts,  and  the  evaporative 
plane  temperature,  Tp. 
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B.1.1  AIR  REGION: 


1»M  h. 
a  a 


ah 

1*M  f  h  +_iAz) 
a  a  dz 


h  (  T  -  T  )  Az  •  1  N  h  Az  •  1 
v  a  s  v  v 


Performing  an  energy  balance  on  CV  I  and  dividing  the  equation  by  Az*l 


yields: 


N  h  +0  =  M  ^  +  h(T  -T  ) 

v  v  ^as  a  dz  a  s 


(B.l) 


Where  ha,  the  total  dry  air  plus  vapor  enthalpy  on  a  per  unit  dry  air  basis  is  given 
by  Flake  [2]  as: 


h  =  CT  +  Y[C  T ,  +  Ah  +  C  ( T - T .  )  ]  (B.2) 

a  a  1  w  dp  vv  dp  J 


The  energy  input  into  the  control  volume,  given  by  the  evaporative  heat  flux  and  the 
supplemental  air  heat  source,  must  equal  the  sensible  heat  gained  plus  the  heat 
convected  into  the  material. 
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B-1.2  DRY  PRODUCT  REGION 


h  (  T  -  T  )  Az  •  1 
v  a  s ' 


N  (  h  +  A(h  )  )  Az 


N  h  Az  •  1 
v  v 


Performing  an  energy  balance  on  control  volume  n,  cancelling  common 
terms  in  both  sides  of  the  equation  and  dividing  both  sides  by  (Az*l)  yields: 

M  ra(h8),  t-t  m 

h(T  -  T  )  =  N  A(h)  +  ^ —  ]  +  X  — - E  (B‘ 

a  s'  v  v7  b  L  8z  J  s  § 

where  A(hv)  is  the  enthalpy  difference  between  the  vapor  at  the  evaporative 
interface  and  the  vapor  at  the  surface.  Equation  B.3  states  that  the  heat  convected 
into  the  material  must  equal  the  change  in  enthalpy  of  the  vapor  being  convected 
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through  the  dry  region  plus  the  sensible  heat  gain  of  the  dry  region  ,  plus  the  heat 
being  conducted  into  the  wet  region. 

Assuming  that  8  is  known  as  a  function  of  X  for  the  material,  (8,z)  can  be 
transformed  to  (8,X)  by  means  of  the  chain  rule. 


as  _  as  ax  dt  -Nv  as 
3z  "  ax  at  3z  “  m  ax 

s 


(B.4) 


Expanding  the  second  term  on  the  RHS  of  equaton  B.3,  the  resulting  term  with  the 
rate  of  change  of  S  with  respect  to  z  is  related  to  the  change  of  enthalpy  of  the  dry 
region  due  to  the  enlargement  of  the  control  volume.  This  term  is  evaluated  at  Tp 
since  the  enlargement  occurs  at  the  evaporative  interface,  which  is  at  Tp. 

To  calculate  the  change  of  enthalpy  of  the  vapor  going  through  the  dry 
region,  it  will  be  assumed  that  the  vapor  leaves  the  material  at  the  mean  temperature 
of  the  dry  region.  Now  the  energy  balance  equation  reads. 


h(T  -T  )  =  N  C  (T  -T  )  + 
v  a  s  v  v  m  p 


M8C  3T  NCT 
s  s _ v  s  p  oo  +  ^ 


T  -T 
s  r 


(B.5) 


Define  the  relative  depth  of  the  wet  product  as 


Performing  an  energy  balance  on  the  control  volume  yields: 


Appling  the  chain  rule,  simplifing  and  using  the  following  expressions: 
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ax 

s  3z 


=  -N 


a;  _  -as  Nv  as 
az  "  baz  ~  bM  ax 
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h  -  h  =  Ah 

V  f  V 
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equation  B.7  becomes: 
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T  -T 


s _ p 
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N  Ah  + 

V  V 


Nvhs  58 

b  ax 


+  M(CC  +xcw) 


aT 

_ p 

dz 


(B.10 

The  second  term  on  the  RHS  of  Equation  B.10,  is  the  enthalpy  associated  with  the 
dry  product  leaving  the  control  volume  because  of  the  evaporative  front  This 
enthalpy  will  be  evaluated  at  Tp. 


B.2  SUMMARY  OF  EQUATIONS  1 1SED  BY  THE  DRYER  CODE 


Energy  Balance,  Air  Region,  CVI 


8h 


Nh+Q  =M^-  +  h(T-T) 

V  V  ^as  a  9z  v  a  s  ’ 


Energy  Balance,  Dry  Region,  CV  H 


M5C8T  NCT 

Q  =  N  C  (T  -T  )+— - -  —01  -  v  s  P  4—  +  Q 

^  ir  u  '  m  n  7  L  ^  L  Av 


v  v  m  P  b  8z 


b  ax 


Energy  Balance,  Wet  Region,  CV  III 


T  • T  Nvh,  38 


aT 


s  p 

K  — „  +  Q„  =  N  Ah  +  .  .  v  .  —  v 

S  5  g  v  V  b  ax  sv  ^  s  W  dz 


+  M„(C c +xc 


Moisture  Mass  Balance,  Air  Region,  and  Product  Region,  CV  I  ,n  and  III. 


N  -M 

v  a  dz  s  dz 


Humid  Air  Enthalpy  Expression 


h  =  CT  +  Y[C  T.  +  Ah  +C  ( T  -  T .  )  ] 

a  a  L  w  dp  v  vv  dp  J 


Dew  Point  Temperature  Definition 


Ysat(Tdp)  =  Ya 


Drying  Rate  Equation  for  Convective  Drying. 


Nv=/rKoDln[^i^.] 


Drying  Rate  Equation  for  Hybrid  Drying. 


% 

Nv=^r^ 


Surface  Energy  Balance  Defining  Wet  Bulb  Conditions 

h(Ta-Tw)  =  Ko  D  In  hfg 

Relationship  Between  Saturated  Humidity  and  Temperature 

Ysat =  Ysat(Ta) 

Characteristic  Drying  Curve  Function 

fr  =  /(<&) 

Evaporative  Plane  Position,  5, for  convective  drying 


-/[-if 


Evaporative  Plane  Position,  6,  for  dielectric  drying 


5_  =  8*  =  ( 1  -  — ) 

b  v  X  ' 
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